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COEFFICIENTS 

HIDEO DEGUCHI* AND MICHAEL OBERGUGGENBERGERt 

Abstract. This article addresses linear hyperbolic partial differential equations with non-smooth 
coefficients and distributional data. Solutions are studied in the framework of Colombeau algebras 
of generalized functions. Its aim is to prove upper and lower bounds for the singular support of 
generalized solutions for wave equations with discontinuous coefficients. New existence results with 
weaker assumptions on the representing families are required and proven. The program is carried 
through for various types of one- and multidimensional wave equations and hyperbolic systems. 
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1. Introduction. This paper is an advanced part of a program that aims at 
studying propagation of singularities and regularity of solutions to linear hyperbolic 
equations and systems with non-smooth coefficients in algebras of generalized func¬ 
tions. To explain the setting, consider a wave equation of the form 

dfu(t, x ) — div (c(f, x) Vu(t, x)) =0, t > 0, x £ R d , 
m (0, x) = uq(x), 9 t w(0, x) = u\(x), x £ R d . 

We are interested in the case of non-smooth coefficients c(t , x) and initial data distri¬ 
butions. The reasons to get into algebras of generalized functions are twofold. First, 
it is well-known that the existence of classical solutions requires a minimal degree 
of smoothness of the coefficients, as does the application of the classical machin¬ 
ery of pseudodifferential operators and Fourier integral operators [5, 6, 22, 40]. In 
intended applications in seismology and continuum mechanics, the coefficients might 
have jumps at irregular interfaces or might be continuous, but non-differentiable paths 
of a stochastic process, thereby lacking the kind of smoothness required in the classi¬ 
cal theory [9, 36]. Second, as the solution is expected to be as singular as the initial 
data, equation (1.1) involves a nonlinear interaction of coefficient singularities with 
singularities in the solution, hence a product of distributions. 

There is a long history of results (e.g. [18, 23, 30, 31, 33, 34, 35]) that equations of 
this type can be uniquely solved in the Colombeau algebra Q of generalized functions. 
This is an algebra containing the space of distributions as a subspace. Non-smooth 
coefficients and distributional initial data can be imbedded in this algebra, and hence 
generalized solutions are available in this setting. Members of the Colombeau algebras 
are equivalence classes of families (w £ ) eg (o,i] of smooth functions that are characterized 
by certain asymptotic properties as e j. 0. Representatives of the Colombeau solution 
satisfy 

dfu e (t, x) — div (c e (t, i)Vu £ (i, x)) = h E (t, x), t > 0, x £ R d , 
u e (0, x) = U 0e (x) + h 0s {x), d t U e ( 0, x) = Ui e (x) + h le (x), X £ R d 
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where h E , ho s , h\ e belong to a certain ideal. The question of regularity and propagation 
of singularities can be cast in the following form: 

(Q) Is it possible to detect the singular support of the generalized solution from 
the asymptotic properties of its representing families (zi E ) eg ( 0 u as e j, 0? 

In this respect, a subalgebra Q°° of Q has been introduced in [34] which satisfies 
17“ fl D' = C°°. It plays the same role in the (/-setting as C°° does in the In¬ 
setting. In fact, the essential steps of regularity theory and propagation of singularities 
for Colombeau solution have been taken in the past decades, when the coefficients 
are regular (C“ or Q°° with some additional properties). This starts with elliptic 
and (/“-hypoclliptic regularity, regularity of solutions with (/“-regular initial data, 
propagation of the (/“-singular support in the case of constant and of (^“-coefficients. 
Further, a calculus of pseudodifferential operators and Fourier integral operators has 
been developed in the Colombeau setting that admits to deduce a number of results on 
the (/“-wave front set of the Colombeau solution when the coefficients are (/“-regular 
[10, 11, 12, 13, 14, 15, 16, 17, 21, 24, 25, 28, 32, 34, 37], 

The central issue of the paper concerns non-smooth coefficients. When imbedded 
in the Colombeau algebra, they should not belong to the subalgebra Q°° and it should 
be possible to detect the singularities that they produce in the generalized solution. 
The difficulty arises that in most of the previous existence results in the Colombeau 
setting, coefficients are required to satisfy certain additional logarithmic or slow scale 
estimates as e j. 0, which makes them (/“-regular. One way out is to change the 
scale on which regularity is described, as was done e.g. in [26]. Another way out is to 
consider only time-dependent coefficients. A number of results on propagation of sin¬ 
gularities for scalar transport equations (including refined properties of the (/“-wave 
front set at characteristic curves with kinks) and on one-dimensional wave equations 
with time dependent coefficients have been obtained in the literature [7, 8, 35]. What 
has been lacking so far is the important case of x-dependent coefficients. 

At interfaces of discontinuity, incoming singularities produce a reflected and a 
transmitted wave. In view of the intended applications in seismology and continuum 
mechanics, it is crucial to be able to detect the (/“-singularities carried by these 
waves. The present paper closes this gap in the literature by determining precise 
upper and lower bounds on the (/“-singular support of the solution in a number of 
paradigmatic cases. 

The paper has essentially two parts. In the first part, existence results on wave 
equations of the type (1.1) in any space dimension are given that do not require a 
logarithmic regularization scale for the coefficients. Similar results are proven for the 
one-dimensional, non-conservative case, one-dimensional hyperbolic systems, together 
with a new existence result (in the algebra (7) for multidimensional wave equations 
with time dependent coefficients. In the second part, upper and lower bounds on the 
(/“-singular support of the generalized solution to the one-dimensional wave equation 
with an initial point singularity and an ^-dependent coefficient with a jump singularity 
are given. The same is done for rotationally symmetric solutions to multidimensional 
wave equations with t-dependent coefficients. The paper is complemented by com¬ 
paring the generalized solution with the distributional solution of the corresponding 
transmission problem, when it exists. The paper shows that in paradigmatic cases the 
Colombeau solution carries the required information on the propagation of singulari¬ 
ties and indeed answers question (Q) affirmatively in these cases. This indicates that 
the Colombeau setting may be appropriate for studying wave propagation phenomena 
in even more singular media, where no classical counterpart exists. 
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The plan of the paper is as follows. In Section 2, required concepts of the theory 
of Colombeau generalized functions are recalled. Section 3 is devoted to new existence 
and uniqueness results for wave equations with non-smooth coefficients in Colombeau 
algebras. The main goal is to show that in many situations bounds on the coefficients 
from above and away from zero suffice - no logarithmic asymptotic estimates are 
needed. In Subsection 3.1, this is done for the conservative wave equations with 
x-dependent propagation speed in any space dimension, in Subsection 3.2 for the 
non-conservative wave equation in one space dimension, in Subsection 3.3 for the 
wave equation with time-dependent propagation speed in any dimension. Subsection 
3.4 on scalar transport equations paves the way for the results on hyperbolic systems 
in one space dimension and x-dependent coefficients in Subsection 3.5. Subsection 3.6 
is a side remark on singularities produced by a discontinuity of the coefficient at time 
t = 0, if compatibility conditions are lacking, and that this classical phenomenon can 
also be detected in the Colombeau setting. 

Section 4 aims at proving upper and lower bounds on the singularities of the 
Colombeau solution to wave equations, given initial data with a point singularity. 
Subsection 4.1 treats the one-dimensional case with a coefficient having a jump at 
x = 0 and with a delta function at x = — 1 as initial data. We prove that the reflected 
and the transmitted singularity can be detected in terms of the t?°°-singular support 
of the Colombeau solution, provided the height of the jump satisfies certain bounds. 
It is also shown that there is no reflected singularity if the coefficient is C/°°-regular. 
Subsection 4.2 addresses the one-dimensional wave equation with a f-dependent co¬ 
efficient having a jump at t = 1. As a preliminary step, a regularity result is proven 
that states that if the initial data are t/°°-regular, then so is the solution off the line 
t = 1. Then we prove an upper bound on the singular support of the solutions for 
arbitrary Colombeau initial data which are t7°°-regular off the origin, generalizing our 
results from [8]. The upper bound consists of the rays emanating from the origin at 
t = 0 and the refracted and transmitted rays starting at the singularity at t = 1. (In 
[8] it was shown that the bound is attained for delta function initial data.) Next, 
we give sufficient conditions on the initial data such that singularities propagate out 
from the origin and remain present along the transmitted rays. Again, this is com¬ 
plemented by precise bounds on the singular support of the solution if the coefficient 
is £/°°-regular. Finally, Subsection 4.3 settles the rotationally symmetric case with 
f-dependent coefficient and delta function initial data in any space dimension. The 
Appendix complements the paper by addressing the question of existence of associ¬ 
ated distributions (distributional limits of the family of regularizations defining the 
Colombeau solution). Extending our results from [18], we show that in one space 
dimension, an x-dependent coefficient with a single jump at the origin and delta func¬ 
tion initial data, the Colombeau solution admits the distributional solution of the 
corresponding classical transmission problem as associated distribution. On the side, 
it is seen that the C°°-singular support of the distributional solution coincides with 
the C/°°-singular support of the Colombeau solution. As a final note, we remark that 
most of our existence results are based on inversions of Colombeau algebras, due to 
the energy methods used in the proofs. (Exceptions are hyperbolic systems in one 
space dimension and the one-dimensional wave equation with t-dependent coefficient, 
in which cases the usual Lj^-based algebras suffice.) 

2. Basic concepts from Colombeau theory. We will employ the special 
Colombeau algebra of generalized functions denoted by Q s in [20] (called the simplified 
Colombeau algebra in [2]). However, here we will simply use the letter Q instead. In 


4 


HIDEO DEGUCHI, MICHAEL OBERGUGGENBERGER 


the context of energy estimates, we will also need L p — L 9 -based Colombeau algebras, 
introduced in [3]. This section serves to recall the definition and properties required 
for our purpose. For more details, see e. g. [3, 20, 34]. 

Given a non-empty open subset Q of M", the space of real valued, infinitely 
differentiable functions on D is denoted by C°°(f2), while C°°(D) refers to the subspace 
of functions all whose derivatives have a continuous extension up to the closure of Q. 

Let C°°(f2)( 0,1 l be the differential algebra of all maps from the interval (0,1] 
into C°°(f2). Thus each element of C 00 (f2)^ 0,1 ] is a family (it e ) ee (g u of real valued 
smooth functions on O. The subalgebra £m{££) is defined by the elements (u £ ) £ g(o,i] 
of C 00 (fl)(°’ 1 ] with the property that, for all K <s £1 and a £ Ng, there exists N > 0 
such that 


sup \d a u E (x)\ = 0(e n ) as e 4- 0. 

xEK 

The ideal AT(fl) is defined by all elements (u e ) £( =m u of C 00 ^)^ 0 ’ 1 ] with the property 
that, for all K m SI, a € N 0 " and N > 0, 

sup \d a u E {x)\ = 0{e N ) as e J. 0. 

x£K 

The algebra C?(f2) of generalized functions is defined by the factor space 

G(tt)=£ M {n)/ 'AA(Q). 

The Colombeau algebra G(£l) on the closure of tt is constructed in a similar way: the 
compact subsets K occurring in the definition are now compact subsets of Q, i.e., 
may reach up to the boundary. Since £m{TL) C £m(Q) and C Af{£l), there is a 

canonical map C?(S2) —^ G(£l). However, this map is not injective, as follows from the 
fact that J\f(£l) H £m(Q) ^ A/"(H). 

Turning to L p — L q - based Colombeau algebras, we assume for convenience that 
the open set Q has the strong local Lipschitz property (cf. [1]), in order to have all 
required Sobolev imbedding and extension theorems at our disposal. Let 1 < p, q < oo 
and to £ Z. Concerning Sobolev spaces, we employ the usual notations: 

H m {£i) = W m ’ 2 (£L), and W°°’ p (n) = nDue 
to the Sobolev imbedding theorems, 

fF°°’ p (n) c IT 00,00 (D) 

for all p , and consequently 

W°°’ p {n) C W°°’ q {n) 

for p < q; further, these spaces are differential algebras. 

Let £ p (Cl) be the differential algebra of all maps from the interval (0,1] into 
W°°’ p ( fi). The subalgebra £m, p {££) comprises the elements (u £ ) £ g(o,i] of £ p (fl) with 
the property that, for all a £ Ng, there exists N > 0 such that 

||<9 a M e ||Lp(n) = 0(e _Ar ) as e j, 0. 

Further, denote by A f p , q (£l) the set of elements (u £ ) £ g(o,i] of £m, p {££) with the 

property that for all a £ Ng and N > 0, 

l|5“u e ||L9(f2) = 0(e N ) ase|0. 
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It can be shown [3] that A f p , q (fl) is a differential ideal in £m, p {G£), thus 

G p ,q(fl) = A/Jj,g(fi) 

is a differential algebra. It has the following properties: (a) For 1 < pi < p 2 < oo and 
1 < q < oo, 


GpugW C S pa ,,(«). 

(b) There exists a canonical mapping from Q p i9 (H) into C?(fl), which, however, is not 
injective. 

Since fl is assumed to have the strong local Lipschitz property, W°° ,p (fl) C 
C°°(fl). It follows that £ p (fl) = £ P (fl), thus replacing fl by fl makes no difference in 
the L p — L 9 -based Colombeau algebras: 

G P , q {fl) = G P: q(fl). 

This holds, in particular, when fl is a strip (0,T) x R d . If fl is bounded, we have in 
addition 


G P ,q(fl) = G P ,q(fl) = G(fl), 

because then all L p -norms are equivalent with the L°°-norm, due to the Sobolev 
imbedding theorem and the boundedness of fl. 

Summarizing, one should keep in mind that G(fl) is defined by local estimates, 
G(fi) by semi-global estimates (up to the boundary, but local on unbounded parts of 
fl), while Q Pt q(£l) is defined by global estimates. 

We use capital letters for elements of G to distinguish generalized functions from 
distributions and denote by (w e ) E6 (o,i] a representative of U € G- Recall that for any 
U,V £ G and a £ Ng, the partial derivative d a U is defined as the class of (9 Q M e ) ee ( 0i i] 
and the product UV is defined as the class of (M e u e )ee(o,i]- Also, for any U = class 
of (u E (t, aO)e€(o,il € (7([0,oo) x R), we can define its restriction U |t=o £ f?(R) to the 
hyperplane line {t = 0} to be the class of (u e (0, a;)) e e(o, 1 ]■ 

Restrictions to open subsets are defined similarly. The support supp U of a gener¬ 
alized function U £ C?(S1) is the complement of the largest open subset of fl on which 
U vanishes. 

Remark 1. For U £ G p , q (fl) and V £ Goo, q {fi), the product UV is a well-defined 
element of G P , q (fl)- 

Remark 2. The algebra G{fl) contains the space £'{fl) of compactly supported 
distributions. In fact, the map 

/ ^ i(/) = class of (/ * p e |n) E g(o,i] 
defines an imbedding of £'{fl) into G(fl), where 



and p is a fixed element of <S(R n ) such that / p(x) dx = 1 and / x a p(x) dx = 0 for 
any a £ Ng, |a| > 1. This can be extended to an imbedding i of the space V'{fl) of 
distributions (the extension is unique as a sheaf morphism). Moreover, this imbedding 
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turns C°°(fl) into a subalgebra of By an anologous imbedding, H/ _00 ’ P (R") is 

contained in Q p<q (R n ), and B roo,9 (IR”) is a subalgebra of Cy P!g (R n ) if p < q. 

Definition 2.1. A generalized function U £ C?(fl) is said to be associated with a 
distribution w £ D' (f2) if it has a representative (u e ) ee (o,il £ £m(£ 1) such that 

u e —^ w in 2?'(fl) as e I 0. 

We write U ~ w and call w the associated distribution of U provided U is associated 
with w. 

Regularity theory for linear equations has been based on the subalgebra G°°(Q) of 
regular generalized functions in CJ(fl) introduced in [34]. It is defined by all elements 
which have a representative (it e )ee(o,i] with the property that, for all K <s fl, there 
exists TV > 0 such that, for all a £ Nq , 

(2.1) sup |<9“'u £ (a;)| = 0(e~ n ) as e f 0. 

xGK 

We observe that all derivatives of u E have locally the same order of growth in e > 
0, unlike elements of £m{ fl)- This subalgebra has the property t7°°(fl) fl 

i(X>'(fl)) = C°°(fl), see [34, Theorem 25.2]. Hence, for the purpose of describing the 
regularity of generalized functions, plays the same role for C/(fl) as C°°(fl) does 

in the setting of distributions. An element U £ G(£l) is said to be ^“-regular on an 
open subset ui of fl, if its restriction U\ u belongs to <5°°(w). An element U £ G P , q { fl) 
is said to be ^“-regular on w, if the canonical image of U in G(w) belongs to G°°(w), 
or - equivalently - if U has a representative (ite) e e(o,i] i n £m,p( fl) which satisfies the 
estimate (2.1) on every compact subset K of u. In either case, the C/°°-singular support 
(denoted by singsupp^oo) of a generalized function is defined as the complement of 
the largest open set on which the generalized function is regular in the above sense. 

We shall also make use of the global versions of regular generalized functions by 
Gff q { fl) defined by having a representative for which there exists TV > 0 such that for 
all a £ Nq, 

II<9“u £ || lp( q) = 0(e~ n ) as £ 4_ 0. 

It is clear that every element of Gff q ( fl) is ^“-regular in the previous sense, but a 
CJ°°-regular element of G P ,q( fl) need not belong to the global space G™ q { fl)- 

We end this section by recalling various specialized asymptotic properties of nets 
and of generalized functions. A net (i , e) e e(o,il is called a slow scale net if 

\r e \ p = 0 (e^ 1 ) as e l 0 
for every p > 0. It is called a log-type net if 

|r £ | p = 0(| log£|) as £4.0. 

Further, if |r £ | p = 0(1) as £ 4- 0, the net is called bounded. We refer to [28] for a 
detailed discussion of slow scale and of logarithmic nets. 

Finally, a generalized function U from G( fl) or from G P ,q{ fl) is called of bounded , 
of logarithmic , or of slow scale type , if it has a representative (Me)eg(o,i] such that 
sup xgn \u E {x)\ is a bounded, log-type, or slow scale net, respectively. If such an 
estimate holds on every compact subset of fl, the generalized function is termed to 
be of locally bounded, locally logarithmic, or locally slow scale type , respectively. 
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Example 1. (a) Let <p be a fixed element of C“(R) such that <p is symmetric, 
<p' > 0 on [—1,0], supp</? C [—1,1] and f R p(x)dx = 1. Put ip e (x) = p{x/e)/e. Then 
U £ Q (R) defined by the class of (<^ e ) ee ( 0 ,i] is associated with the delta function, and 
sing supple*, U = {0}. 

(b) On the other hand, if U £ (/(R) is defined as the class of (<p?i(e))ee(o,i], where 
h(e) > 0 and (l//i(e)) E g(o,i] is a slow scale net, then it is associated with the delta 
function again, but singsupp^oo U = 0. In particular, this holds for the logarithmic 
regularization h(e) = 1/| loge|. 

At this stage, we are in the position to clarify the relations of the singular support 
of a Colombeau generalized function with the classical singular support of associated 
or imbedded distributions. First, the relation 0 l{V M) = c°° (w) holds for 

every open subset u of f2, whence 

sing suppgoo i(f) = singsupp / 

for every distribution / £ V'(Q). However, if U £ £/(fl) admits w £ as associ¬ 

ated distribution, there is no general relation between the (/“-singular support of U 
and the singular support of w. Indeed, the generalized function U from Example 1(b) 
provides an example of a generalized function whose (/“-singular support is empty, 
while its associated distribution has nonempty singular support. On the other hand, 
the (/“-singular support of the generalized function V represented by (£¥’e( a; ))Ee(o,i] 
is {0}, while its associated distribution (which equals zero) has empty singular sup¬ 
port. Consequently, knowledge of the singular support of the associated distribution 
does not suffice to determine the (/“-singular support of a Colombeau generalized 
function. 

Following [32, 37], an element U £ 17(H) is said to be strongly associated with 
a distribution w £ X>'(f2), if for every ip £ X>(fi), |/ u E (x)ip(x)dx — (w,ip)\ = 0(e p ) 
for some p > 0. Further, U and w are said to be equal in the sense of generalized 
distributions , if this estimate holds for every p > 0 (see [4]). If U is strongly associated 
with w, then singsupp w C singsupp^oo U. If U and ui are equal in the sense of gener¬ 
alized distributions, then singsupp ui = singsuppgoo U. However, solutions of partial 
differential equations are rarely known to be strongly associated with a distribution. 

As a last ingredient, we need the notion of a compactly bounded or c-bounded 
generalized function (see [20, Def. 1.2.7]). Only the scalar case is required here. An 
element U £ £7(R) is called c-bounded, if it has a representative (u e ) E e(o,i] with the 
property 

(2.2) \/K <s R 3K' <£ R such that u e {K) C K' 

for all sufficiently small e. The importance of the notion lies in the fact that compo¬ 
sition of c-bounded generalized functions is possible. More precisely, if V £ (?(R) and 
U is c-bounded, then V o U is a well defined element of t?(R). 

3. Existence and uniqueness of generalized solutions. In order to be able 
to detect the (/“-singular support of a generalized solution to the wave equation 
with non-smooth coefficients, it is necessary to have existence and uniqueness results 
that do not require logarithmic bounds on the coefficients (as is the case with the 
classical existence results for hyperbolic systems in Q). Indeed, if the (discontinuous) 
coefficients of the equation are regularized by convolution with a modifier (</ 5 /i(£))ee(o,i] 
as in Example 1(b), they become (/“-regular generalized functions and contribute no 
(/“-singularities to the solution, in particular, no reflected or refracted waves, as will 
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be seen below. It is the purpose of this section to provide existence and uniqueness 
results that can dispense with such slow scale or logarithmic regularizations and get 
along with standard modifiers (</j E ) Eg (o,i]■ 

Indeed, there are a number of cases in which no logarithmic estimates are needed, 
when the setting of C/ 2,2 is used: the conservative wave equation with propagation 
speed depending on x only (in any space dimension), the one-dimensional wave equa¬ 
tion with propagation speed depending on x , and the wave equation in any space 
dimension with t-depending propagation speed, bounded and of bounded variation. 
(The one-dimensional, t-dependent case was settled in [8].) In all cases, a bound away 
from zero is required as well. 

An alternative approach would consist in measuring the regularity of a gener¬ 
alized function based on an /i(e)-scale, replacing the d(e: _Ar )-bound in (2.1) by an 
0(/i(e) _Ar )-bound, as has been suggested in [26]. This requires choosing h(e) appro¬ 
priately for the specific equation at hand. In this paper, we aim at removing the need 
for introducing a scale other than the standard scale h{e) = e in the existence results, 
in order to be able to apply the standard t/°°-regularity concepts. 

3.1. Existence theorems for conservative wave equations. We start with 
the conservative wave equation 

dfu{t, x) — div (c(x)Vw(t, a;)) =0, t > 0, x £ K d , 
u(0, x) = uq(x), dtu(0, x) = u\(x), xGM. d . 


When considering solutions whose representatives satisfy L 2 -estimates in t and x, we 
shall have to work on strips [0, T] x R d , but with arbitrary T > 0. Thus looking for 
generalized solutions in the Colombeau algebra £ 2 , 2 ([0, T] xl d ), we interpret Equation 
(3.1) as 

dfU — div (CS7U) = 0 in £2,2([0, T] x R d ) 

U\t=o = U 0 , d t U\t=o = Ux in £ 2 , 2 (M d ) 

where the coefficient C is a generalized function depending on x £ only. We make 
the following assumptions about C. 

(Al) C belongs to f?oo, 2 (® d )> that is, it has a representative (c E ) eg ( 0 ,i] all whose 
derivatives have global moderate L°°-bounds; 

(A2) C is strictly positive, that is, there is m > 0 such that, for every representa¬ 
tive, 0 < e m < inf xgR d c e (x) as e i 0. 

Theorem 3.1. Assume that Uq,Ui £ I/ 2 , 2 (R d ) and that the assumptions (Al) 
and (A2) hold. Let T > 0. Then the conservative wave equation (3.2) has a unique 
solution U £ f/ 2 , 2 ([0,T] x R d ). 

Proof. Given representatives of Uq , U\ and C, classical existence theory provides 
a prospective representative (n e ) E g(o,i] where each u e belongs to H°°([0,T] x K d ). It 
satisfies the equation 

dfu e (t,x) - div (c e {x)Vu e {t,x)) = 0, t £ [0,T], x £ R d , 

U e (0, x) = UQe(x), d t u e (0,x) = u le (x), X £ R d . 


We have to prove that (u £ ) eg ( 0l i] and all its derivatives are moderate in the L 2 -norms. 

Multiply the differential equation (3.3) by dtu e and integrate by parts (with re¬ 
spect to x): 



u e d t u e + c e (*)V 


u 


£ 


■ V<9 t u E ) dx = 0, 
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Id 
2 dt 


f (\d t u E \ 2 + c E (x)\Vu e \ 2 ) dx = 0. 
J R d 


This gives a moderate bound on the L 2 -norms of dtu e and Vit e in terms of the initial 
data. Since u £ (t, x) = itoe + f 0 dtu e (s, x) ds , this also entails moderate bounds on the 
i 2 -norm of u E on the strip [0, T} x K d . 

Next, dtu e satisfies the wave equation as well, thus the same argument gives 
moderate bounds on the L 2 -norms of d 2 u e and 9 t Vu e . From the equation we obtain 


(3.4) 


dfU E — c f Au f — Vc F • Vu F = 0 


and in turn moderate bounds on the L 2 -norm of Ait e , involving also the L°°-norm of 
Vc £ . But ||M e || L 2 ( R d) + ||Vu e || L 2 (R d) + ||Aw e || i 2 (R d) is an equivalent norm on H 2 (R d ), 
thus a moderate bound on the H 2 -norm of u E on [0,T] x is obtained. 

Next, we are aiming at moderate bounds on the F/ 3 -norm of u E . Considering that 
d 2 u E satisfies again the wave equation (3.3), one gets moderate bounds on dfu E and 
S7d 2 u E . Denoting by d one of the x-derivatives d Xi , we get from (3.4): 

d 2 du e — c E Adu E — (dc e )Au e — V9c e ■ Vu e — Vc e • Vdu e = 0. 

All terms except c e A du e are already known to have moderate bounds by the previous 
steps, thus the equation entails a moderate bound on the L 2 -norm of A du e , and in 
turn on the H 2 -norm of du E . It remains to estimate terms of the form d Xi d Xj d t u e . 
But an estimate of A d t u e follows again from (3.4), and the equivalent-norms-argument 
completes the step. Collecting all terms, we see that the H 3 -norm of u E has moderate 
bounds. In the same way, the higher order derivatives can be estimated. 

To prove uniqueness, we consider the inhomogeneous equation 


d 2 u e {t, x) — div (c e (x)Vu e (t, x)) = h e (t, x), f€[0, T], x£ 

it E (0,x) = ho e (x), d t u e (0,x) = /n e (x), x e M d 


(3.5) 


where h E belongs to A/ 2 , 2 ([0,X 1 ] x ffir) and ho E: hi e belong to A/ 2 , 2 (K“)- We have to 
prove that u E belongs to A/^QO, T] x R d ) as well. Thanks to [10, Prop. 3.4] it suffices 
to prove the negligibility of the L-norm of u E (no derivatives needed). As before, we 
get 


1 d 

2 dt 


'r^ 


(\d t u E \ 2 + c E (x)\Vu E \ 2 ) dx = / h E d t u E dx. 

JR d 

Integrating, inserting the initial data and using the nonnegativity of c E gives 

\\dt u e(t, •)lli 2 (Rrf) A II|| L 2 (R d ) + II c e || (]R d ) II II z. 2 (]R d ) 

+ / \\he{s,-)\\‘L 2 (m. d )d’S + / |j^M £ (s,-)|| 2 2( R d) ds. 
Jo Jo 


By assumption, all the terms on the right hand side except the last one are negligible. 
By Gronwall’s inequality, the negligibility of |||| z, 2 ([ 0 ,^] xK d ) an d then in turn of 
IMlL 2 ([ 0 ,T]xRrf) follows. 0 

If the initial data are ^“-regular, so is the solution, provided the coefficient 
C £ £ 00 , 2 (R d ) satisfies slow scale estimates, more precisely, assume that: 

(A3) C and all its derivatives are of slow scale type, that is, there is a representative 
(c e ) e6 (o,i] all whose derivatives have global L°°-bounds in terms of a slow scale 
net; 


10 


HIDEO DEGUCHI, MICHAEL OBERGUGGENBERGER 


(A4) C is slow scale invertible, that is, there is a slow scale net (r e ) e6 (o,il such 
that, for every representative, 0 < l/r £ < inf^ggd c £ (x) as e 4- 0. 

Corollary 3.2. Assume that C satisfies the assumptions (A3) and (A4). Let 
Uq, U\ £ Then the solution U of (3.2) belongs to ^ 2 ^ 2 ([0, T] x M* 1 ). 

Proof. It is quite clear by following the same arguments. At each step, products 
of derivatives of c e and u e occur (when computing the initial data and when solving 
for the highest order derivative), containing only two factors each. Take pj > 0, 
YlpL 1 P.i = 1; where each pj corresponds to one of the recursion steps. Any derivative 
of c £ is of order 0(e~ Pj ) for whatever j. At each step, the previously estimated terms 
are multiplied by a factor of order 0(e~ Pj ). Hence all derivatives of u E will be of order 
0(e _Ar_1 ), where 0{e~ N ) is the order of the ^“-initial data. 0 
Remark 3. The equation 

dfu(t,x) — div (c(x)Vw(i, x)) + $(x)u(f, x) = 0, t £ [0, T], x £ R d , 

^ ^ u(0,x) = uo(x), dtu(0, x) = iti(x), xeR d 


with a non-negative potential $ belonging to { 7 oo, 2 (R d ) can be treated in the same 
way. 

Remark 4. In the non-conservative case 
^ ^ d?u E (t,x) - c e (x)Au e (t,x) = 0, t £ [0,T], x £ R d , 

M £ (0, x) = Uq s (x), d t u £ (0, x) = Uie(x), X £ K d 


the energy estimate leads to 


and so 

(3.8) 



Rt'Ue T Cg(x) Tl-l £ 


Vd t u e + Vc £ • Vu e d t u e ) dx = 0 


1 d 

2 dt 


(\d t u E \ 2 + c e (x)|Vu £ | 2 ) dx < f |Vc e | (|Vw e | 2 + |9 t u e | 2 ) dx. 


If the initial data and the L°°-norm of |Vc £ | are of order e N , an application of the 
Gronwall inequality results in a bound on the ffCnorm of u e of the form e~ N e Te , 
which clearly is not moderate. We see that logarithmic estimates on Vc £ would be 
needed to get a moderate bound on the .ffCnorm of u E , at least when continuing with 
this energy estimate. 

This problem can be circumvented in the case of one space dimension. 

3.2. Existence theorems for non-conservative wave equations. We con¬ 
sider the one-dimensional non-conservative wave equation 

(3.9) = t >0,* £R , 

u(0, x) = «o(x), dtu(f),x) = u\{x), x £ R. 


Taking one derivative with respect to x results in the conservative wave equation 

dtd x u{t,x) - d x (c{x)d x d x u{t,x)) =0, t > 0, x £ R, 
m(0,x) = Uq(x), d t d x u(0,x) = u[(x), x£K 

for d x u. Thus the results from the previous subsection can be applied, and this will 
be the path around the difficulty mentioned in Remark 4. Note that this argument 
relies on the fact that the space dimension equals one. 
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Again, we interpret Equation (3.9) as 


(3.11) 


d/U - Cd 2 x U = 0 in g 2t2 ([ 0, T] x R) 

U\t=o = U 0 , d t U\ t= o = Ui in fy 2 , 2 (M) 


where the coefficient C is a generalized function depending on x £ R only. The 
assumptions on the coefficient C £ f?oo, 2 (R) are the same as in Subsection 3.1 with 
71=1. 

Theorem 3.3. Assume that Uo, U\ £ 1?2,2(R) and that the assumptions (Al) and 
( A2 ) hold. Let T > 0. Then the non-conservative wave equation (3.11) has a unique 
solution U £ 02,2([0,T\ x R). 

Proof. We begin by constructing a prospective representative (w e ) e g(o,i] of the 
solution U. Observe that the classical solution u e to (3.7) belongs to H°°([ 0, T] xR), as 
follows from applying the energy estimate from Remark 4 to all derivatives. As lined 
out in Remark 4, we cannot use the energy estimate (3.8) directly to obtain moderate 
bounds on the iJUnorm of u E . Rather, we employ the observation that d x u E (t,x) 
solves the conservative wave equation (3.10). As shown in the proof of Theorem 3.1, 
(9a;We)ee(o,i] belongs to £m,2([0 ,T\ x R). It follows from the wave equation (3.7) and 
assumption (Al) that 


d?u e (t,x) = c e (x)d x u e (t,x) 

belongs to £m,2([0,T] x R) as well. Taking into account the initial data 
d t u E (0, x) = u le {x), u e (0, x) = u 0e {x), 


a double integration finally shows that (rt £ ) ee (o,il belongs to £m,2([0,T] x R), and 
hence may serve as a representative of the generalized solution U. 

Uniqueness is proven along the same lines. □ Global regularity follows immediately 
from Corollary 3.2, because a generalized function U has the £7°°-property, whenever 
its derivatives from a certain order onwards have the C?°°-property. 

Corollary 3.4. Assume that C satisfies the assumptions (A3) and (A4). Let 
Uq, U± £ £? 2 ° 2 (®)- Then the solution U of { 3.11) belongs to f7“ 2 ([0, X 1 ] x R). 

Remark 5. Construction of solutions in f/([0, T] x R d ). Every global £m, 2 -estimate 
implies local fM-estimates. Thus the results of the previous subsections show that 
Equations (3.2) and (3.11) possess solutions in f?([0, T] xR d ), respectively in <?([0, T] x 
R), provided the initial data have a representative in £m, 2 - This is the case, for 
example, if the initial data have compact support. 

However, uniqueness is lost when fy 2 ,2 is replaced by Q. As an example, consider 
the wave equation 



u E (Q,x) =(p{x-1/e), d t u e { 0,x) = 0 


where <p £ T>{ R) with y>(0) = 2, say. This net of functions represents the zero element 
in £?(R). Nevertheless, the corresponding solution 


u e (t,x) = ~((p(x - 1/e - t/e) + <p(x - 1/e + t/e)) 


does not belong to A/"([0,T] x R). Indeed, it e (l,0) = 1 for all sufficiently small e. 
Note, however, that ip(x — 1/e) does not belong to A/ 2,2 (®)j hence does not represent 
the zero element of ^ 2 , 2 (®), for which case the solution has been proven to be unique. 
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3.3. The time-dependent case. This subsection is devoted to the wave equa¬ 
tion with coefficient c(t) depending on time only. We consider 

dfu(t,x) — c(t)Au(t,x) = 0, t> 0, i £ I J , 
u(0, x) = uq(x), dtu(0, x) = u\(x), let 11 . 


Looking for generalized solutions in the Colombeau algebra C/ 2j2 ([0, T\ x we again 
interpret Equation (3.12) as 

d 2 U - CAU = 0 in Q 2 2 ([0, T] x R d ) 

(3 13) ’ 1 J 

U\ t= o = Uo, a t U\t=o = U 1 inS 2 , 2 (ffi d ) 


where the coefficient C is a generalized function depending on t £ [0, 00 ) only, about 
which we make the following assumptions. 

(A5) C belongs to <?oo, 2 [0, 00 ), that is, it has a representative (c £ ) £g ( 0il ] all whose 
derivatives have global moderate L°°-bounds; 

(A6) C is bounded away from zero, that is, there is b > 0 such that 0 < b < 
infte[o,oo) c e {t) as e | 0; 

(A7) C' has global id-bounds in the sense that for all T > 0 and some representa¬ 
tive, ||c(.||ii(o,T) is bounded independently of e as e \. 0. 

Theorem 3.5. Assume that Uq,U\ £ f? 2 . 2 (K. d ) and that the assumptions (A5) — 
(A7) hold. Let T > 0. Then the wave equation (3.13) has a unique solution U £ 
£ 2 , 2 ([0,T] xR d ). 

Proof. Given representatives of do, U\ and C, classical existence theory provides 
a representative (u £ ) £g (o,i] where each u e belongs to H°°([0,T] x R d ). It satisfies the 
equation 

14 ^ 9fU e (t, x) — c s (t)Au E (t, x) = 0, t £ [0,T], x £ R d , 

u e { 0, x) = u 0e (x), d t u E ( 0, x) = ui E (x), x £ R d . 


In order to prove that (it £ ) £g (o,i] and all its derivatives are moderate in the I 2 -norms, 
multiply the differential equation (3.14) by dtu e and integrate by parts (with respect 
to x): 


/ (d?u E d t u E + c e (t)\7u e ■ S7d t u E ) dx = 0. 
J R d 


Observing that 

d t (c E (t)\Vu e \ 2 ) = c E (t)d t \Vu E \ 2 + c' e {t)\Vu E \ 2 = 2 c E (t)Vu E ■ S7d t u e + c' e (t)\X7u E 
we arrive at 

~ [ (\dt u e\ 2 + c E (t)\Vu E \ 2 ) dx = c' e (t)^- j \Vu e \ 2 dx 
z ai j K d z j R d 


whence 


/ R d 


\Vu E \ 2 dx < \ [ (|ui e | 2 +c £ (0)|Vu Oe | 2 )dx-l- [ [ \Vu e \ 2 dxds. 

0 Jm. d Jo b Js. d 


By Gronwall’s inequality, 

[ \Vu E \ 2 dx <j I {\u lE \ 2 +c £ (0)|Vu 0£ | 2 )dxexp( f 
J K d 0 J R d Jo 




ds ) 
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which is finite by assumption. From 


\d t u E \ 2 dx < / (\ui E \ 2 + c E (0)\Vu OE \ 2 )dx + |c'(s)| / \Vu e \ 2 dxds 


we get a bound on the U 2 -norm of dtu e , and then in turn on u e on the strip [0, T\ x R d . 

Next, observe that all ^-derivatives of u E satisfy the same equation (3.14), so sim¬ 
ilar bounds are obtained on their £ 2 -norms. Bounds on the t- and mixed derivatives 
can be obtained successively from the equation. 

Uniqueness is proven in the same way. □ 

Remark 6 . (i) As can be seen from the proof, it would suffice to have a lower 
bound on c E in terms of the inverse of a logarithmic net (in place of the constant b). 

(ii) Usually, the wave equation 

d 2 u(t , x) — c(t) 2 Au(t, x) = 0, t > 0, x € R d 


is written with c(t) 2 in place of c(t). In this case one has to assume a uniform, global 
upper bound B on c E (t), in addition, in order that the derivative of c E (t) 2 remains 
integrable. 

3.4. Scalar transport equations. Existence and uniqueness results for hyper¬ 
bolic first order systems of partial differential and pseudodifferential equations have 
been obtained in the Colombeau setting in [18, 31, 33, 35], requiring certain logarith¬ 
mic or slow scale estimates on the propagation speed. In the next subsection, we will 
show that these estimates can be replaced by requiring upper and lower bounds. The 
present subsection serves to prepare the arguments in the scalar case. 

We consider the scalar transport equation 


(d t + c{x)d x )u{t,x) = 0 , 
it(0, x) = uq(x) 


with transport velocity depending on x only. We are interested in existence and 
uniqueness of a solution in the Colombeau algebra (7([0, oo) x R). Thus we interpret 
problem (3.15) as 

d t U + Cd x U = 0 in f7([0, oo) x R) 

(3.16) 

U | t=0 = Uq, in G(R) 


where the coefficient C is a generalized function depending on x G R only. In the 
quoted earlier results, the coefficient C was allowed to depend on x and t, but it 
was required that C was of bounded type and d x C of logarithmic type. Here, the 
logarithmic type will be replaced by an estimate from below. We make the following 
assumptions. 

(A 8 ) C belongs to £ 00 , 2 0&), that is, it has a representative (c e ) e€ (o,i] all whose 
derivatives have global moderate I/°°-bounds; 

(A9) C is bounded away from zero and is of bounded type, that is, there are 
constants &o,^i >0 such that, for every representative, 0 < 60 < c E (x) < b\ 
for all x £ R as e 4 - 0 . 

Note that at fixed e the coefficient c E (x) does not change sign on R. The standard 
situation, arising from regularization of piecewise constant functions which are either 
positive everywhere or negative everywhere, is that c E (x ) has the same sign for all x 
and all e. 
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For x G R, define 


C e (x) 



dy 

c E {y)' 


For the sake of exposition, consider the case when c E > 0. Then the inequalities 
0 < bo < c E (x) < b\ imply that 


f <C e (x)<^(x>0), 
o 1 0 0 

j- < O e (x) < ^-(x < 0). 

bo b i 

Thus C e is c-bounded (cf. (2.2)) and surjective. In the case under discussion, it is 
also strictly increasing. Its inverse satisfies 

boy < C~ 1 (y) < biy(y > 0), 
hy < C~ l (y) < b 0 y{y < 0). 

Thus C~ 1 (y) is c-bounded, and its derivative is given by ( C E 1 )\y) = c E [C~ 1 (y)). 
Due to the c-boundedness of C“ 1 (y), this is a well-defined generalized function. 
In conclusion, (C' £ ) e6 (o ! i] represents an invertible generalized function, which is c- 
bounded, together with its inverse. 

The characteristic curves x(r) = 'y E {t,x,r) are obtained as solutions to the dif¬ 
ferential equation 


x{t) =c e {x(r)), x{t) = %(t,x,t) = x 


where the dot signifies differentiation with respect to r. Integrating this equality gives 


or 


x(s) ds 
c e (x(s)) 


= r-t , 


C e (x(t )) -C e {x) =T-t , 

7 E {t,x,r) = x(t) = C~ x (C E {x) + t - t). 

Due to the c-boundedness of all participating generalized functions, this represents 
again a well-defined generalized function, which is c-bounded. 

Theorem 3.6. Assume that Uq € C?(M) and that the assumptions (A8) and (^49) 
hold. Then problem (3.16) has a unique solution U G t?([0,oo) x R). 

Proof. A representative of the prospective solution is u e (t,x) = Uo e ("/ E (t,x,0)). 
Since is c-bounded, u E represents a well-defined element of <7([0, oo) x R). A simple 
calculation shows that 

d x 7 e (f,x, 0 ) = c e ( 7 E (f,x, 0 ))/c e (x), dt^ e (t, x, 0) = -c e ( 7 e (t, x, 0)), 
thus u E solves the differential equation. 

To prove uniqueness, it suffices again to establish the null estimate of order zero. 
Thus assume that some moderate v E satisfies 

(d t + c s {x)d x )v e {t,x) = h e (t,x), 
v E (0,x) = h 0e (x) 
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where h E and hg E represent null elements. It holds that 

v e (t,x) = h Oe (7e(t,x,0)) + / h E (r, %(t, x, t)) dr. 

Jo 

Due to the c-boundedness of 7 e (t, x, t), null estimates of h e and ho e on compact sets 
translate into similar estimates on v s , which hence is negligible. 0 

Remark 7. In the same way, existence and uniqueness of a solution in Sea ,00 ([0, 00 ) x 
R) can be proven. 

Example 2. The boundedness of the coefficient c E (x) does not suffice to obtain a 
moderate representative of the solution, as shown by the following example in which 
c E (x ) changes sign. Indeed, consider equation (3.15) with 


c e (x) = — tanh —. 

e 


Again, the characteristic curves x(r) = j e (t,x,r) are obtained as solutions to the 
differential equation 


x(r) = — tanh 

Integrating this equality gives 
f T x(s) ds 


u{t) 


s{t) = 7 e(t,X,t) = X 


= £ 


t tanh^i£l J x/e tanh y 


:(r)/e 


dy 




or 


log 


/ sinh £M\ 

) 


t T 

£ £ 


and finally 


7 e (f,x,r) = x(t) = eArsinh(e^ T ^ e sinh—). 

£ 

For the sake of simplicity, assume the initial data are classical, i.e., uq £ H°°( R). 
Then the corresponding solution is given by 

(3.17) u e (t,x ) = ito( 7 e(£,x,0)) = Uo(^£ Arsinh (e t//e sinh—)^. 

Observe that 

d x u E (t,x) = u' 0 ^£ Arsinh (e t/,£ sinh— ^1 + (e^ £ sinh—) 2 ^ e t/,e cosh—. 

Thus d x u E (t, 0) = Uo(0)e*/ e is not moderate if Up(0) ^ 0. 

To further understand the behavior of the solutions u E (t,x), observe that they 
have a distributional limit, namely 


u £ (t, x) —»• 


uq(x + t), x > 0 , 
uq(x — t), x < 0 . 
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Indeed, it follows from l’Hopital’s rule that 


lime Arsinh (e tA sinh —) = x + t sign a; 
eJ.0 v e 


for t > 0, x ^ 0. As u o was supposed to be a classical bounded function, the assertion 
follows. 

Example 3. Consider again equation (3.15), taking the coefficient 


c e (x) = tanh ■ 


of the opposite sign. In this case, the solution 

u E (t,x) = uq (e Arsinh (e~ tA sinh— )) 

is moderate, more precisely, defines an element of £m([0,oo) x R). This can be 
most easily seen by observing that all partial derivatives of the function v(t, x ) = 
Arsinh (e _t sinhx) are bounded and that u £ (t,x) = uo(ev(t/e, x/s)). The solutions 
u e (t,x) again have a distributional limit, namely 


{ uo{x — t), x > t, 

-uo(O), \x\ < t, 

uq(x + f), x < —t. 

Indeed, for |x| < t, we use that | sinh || < so e Arsinh (e _< ' £ sinh|) A 0 as 

e 4- 0. For |x| > t, Arsinh (e -t / £ sinh |) —> ±oo and l’Hopital’s rule gives again the 
result. 

3.5. Hyperbolic systems. The results on transport equations can be extended 
to hyperbolic (m x m)-systems in diagonal form 

(d t + Ci(x)d x )u,i(t,x) = Y^jLi a ij(t,x) u j(t,x), 

(3.18) 

Ui( 0,x) = u 0i (x) 

for i = 1,... ,m. Concerning the coefficients Cj(x), the same assumptions as in the 
previous subsection are in effect. Concerning the coefficients a,, (t.x). logarithmic 
type has to be assumed as in earlier papers. The novelty of the result is that no 
logarithmic estimates on the derivatives of Cj(x) are required. 

Again, we interpret problem (3.18) in the Colombeau setting as 


dtUi + Cid x U = J2T-1 A ijUj in £([0, oo) x R) 

3.19 3 

Ui\t=o = Uoi, in f?(R) 

for i = 1 ,...,m, where the coefficients Ci are generalized functions depending on 
xel only. 

Theorem 3.7. Assume that the initial data Uoi belong to f/(R), the coefficients Ci 
satisfy assumptions (A8) and (A9), and the coefficients Aij belong to G([0, oo) x R) 
and are of local logarithmic type. Then problem (3.19) has a unique solution u £ 
<7([0, oo) x R). 

Proof. Given representatives of C/o*, Ci and Aij. classical existence theory pro¬ 
vides a representative (n £ ) e g(o.i] of the prospective solution where each u e belongs to 
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C°°([0, oo) x R) and solves equation (3.18), with uoi E , Cj E and aij e in place of u^, Ci, 
a,ij. Integrating along characteristic curves, we may rewrite (3.18) as the system of 
integral equations 



Fix T > 0 and a compact interval K C R and let Kt be the region bounded by 
t = 0, t = T and by the straight lines with slope supCii E (left) and infers (right), 
determined according to condition (A9). If a point (t,x) belongs to Kt, all backward 
characteristic curves of the system connect it to a point in K. Thus we can take 
the L°°-norms on both sides of the system of integral equations on Kt- Gronwall’s 
inequality and the logarithmic type of dij E show that the L°°-norm of Ui E on Kt is 
moderate, i = 1 ,..., to . 

Next, dtUi e satisfies the equation 


T Ciei.x'jdx^dtUieit, x) — — i 0 ’ij E (t, x')df J Uj E (t, x) 

+ YJj=\ d taije(t,x)Uje{t,x), 
dtUi £ (0,x) ^2j=l Q-ijei®, (*e) ^i E (^)'lig^ e: (x). 


(3.20) 


The initial data are still moderate; the second summand in the differential equation 
is known to be moderate from the previous step, and so the Gronwall argument can 
be applied to dtUi e to obtain its moderateness. Finally, equation (3.18) and the lower 
bound on Cj E (a;), due to assumption (A9), entail the moderateness of d x iii e . 

Differentiating the equation again, the moderateness of is obtained in the 

same way. The previous equation (3.20) now entails the moderateness of d x dtUi E . 
Next, we differentiate the system (3.18) with respect to x to obtain 


(3.21) 


T ^is^X^dx^dxUieit, x') — ^^/j= 1 ^ije (^; x')d x 1 lj E (t, x'j C^ e (^x)d x lli E {t,X^ 

T y.j -1 d x dg E (t, x)u-i E (t , x), 
d x u ie ( 0,x) = u' 0iE (x). 

At this stage, all terms in the equation (3.21) have already been proven to be moderate, 
except Ci E (x)d x Ui E (t 1 x), which hence must be moderate as well. Using the lower bound 
on Ci E (x) we arrive at the moderateness of d x Ui e (t, x). 

Next one proceeds step by step: first d^Ui e is estimated by differentiating the 
system of equations once more and using the Gronwall argument, then d^d x Ui E by 
collecting terms in (3.20), differentiated with respect to t, then dtd x Ui E is shown to 
be moderate by collecting terms in (3.21), differentiated with respect to t , and finally 
d x Ui E by differentiating (3.21) with respect to x, and so on. 

Uniqueness is proven along the same lines. □ 

Note that the argument depends essentially on the fact that the coefficients Ci 
depend on one variable only. If the Ci depended on both variables, the recursive pro¬ 
cedure could not be initiated without logarithmic assumptions on the first derivatives 


of Ci. 


3.6. Singularities due to non-compatibility. If the coefficients of a hyper¬ 
bolic equation have a jump discontinuity at some point Xq, a singularity emanating 
from (0,a;o) will be produced, in general. This happens even for smooth initial data, 
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unless suitable compatibility conditions at the point Xq are satisfied. This section 
serves to exhibit this phenomenon in a simple example and to emphasize that it may 
arise in general. In the subsequent sections we will not further elaborate on this 
effect and usually consider only data which actually vanish near the jumps of the 
coefficients. 

For the purpose of illustration, consider the equation 

^ 3 22 ) (d t +c(x)d x )u(t,x) = 0, 

u{ 0, x) = u o(x) 


where c(x ) is piecewise constant with a jump at x = 0. The purpose of this subsection 
is to study the singularity emanating from the origin, both in the classical case and in 
the Colombeau framework, and point out compatibility conditions of the jump data 
with the initial data. In particular, we want to demonstrate in a simple example that 
the Colombeau solution is not ^“-regular along the characteristic curve emanating 
from the origin, in general. 

For the sake of exposition, assume that 


c{x) 


1, x < 0, 

2, x > 0. 



Fig. 1. Geometry of the transmission problem. 

We first recall the classical arguments, interpreting problem classically as a trans¬ 
mission problem. The standard transmission condition is that u(t, x) should be con¬ 
tinuous across x = 0. Observe that in Figurel the backward characteristic through 
(t, x) hits the t-axis at time t — x/2. Thus the solution of the transmission problem is 

{ uo{x — t) in (I), 
uo (—t + x/2) in (II), 
uq(x — 2t) in (III). 

Regions (II) and (III) are separated by the line x = 2 t. It is seen that the solution is 
continuous across this line. Next, 

{ u' 0 {x-t) in (I), 

^u' 0 (-t + x/2) in (II), 
u' 0 (x — 2 1) in (III). 

Thus d x u[t,x) is continuous across x = 2t if and only if Uq( 0) = 0. It is never 
continuous across x = 0, unless u' 0 vanishes on the negative half-axis. 
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Next, we want to show that the same phenomenon is observed in the Colombeau 
setting. We regularize the coefficient c(x) by means of a nonnegative, standard modi¬ 
fier <p e as in Example 1, i.e., ip is smooth, nonnegative, has integral one, is symmetric, 
and its support is contained in [—1,1], p E (x) = ip(x/e)/£. Let C be the class of 
(c * v? E ) eg (o,i] in Given initial data Uq £ 5(K), problem (3.22) has a unique 

solution U £ <?([0,oo) x R). Assume the initial data belong to C/°°(]R). We want to 
show that the singularities across x = 0 and x = 21 can be detected, i.e., that U is 
not Q°° across these lines, in general. 

For the sake of the argument, we take Uo as the C°°-function uq(x) = x and 

/ ;x r x / £ 

Ve{y)dy = 1 + / <p(y) dy. 

-oo J — oo 

Recall from the Subsection 3.4 that the characteristic curves are given by 
7 e{t,x,r) =x{t) =Cj 1 (C e (x)+T-t). 

Due to the simple initial data we use, the solution is given by 

U e {t,X ) = 7 e (t,X,0). 


Recall further that 

<9 x 7 £ (t,x,0) = a t7e (i,x, 0) = —Ce( 7 e (i, x, 0)). 

c e ( x) 

We compute 


dhe(t,x, 0) = - 


c'e( x )Cehe(t,X,0)) c' e (%(t, X, 0))d x %(t, X, 0) 


Ce(xy 




For t > 0 and x near 0, we have that 

c E ('Y e (t,x,0)) = 1, c'(7 £ (t,a:,0)) = 0. 


Thus for those points we obtain 


d x ^e{t,x,G) 


1 

Ce(x)’ 


8he{t,X,0) 


Cs( x ) 

Ce(x ) 2 ’ 


dhe(t,X,0) 


2 «(x)) 2 

c s {x) 3 


c e ( x) 2 ' 


Due to the symmetry assumption on ip we have that c £ (0) = §. Further, c' e {x) = 
c"(x) = |), and so on. Thus if we take <p(0) = a / 0 and p'{0) = 

<p"(0) = ■ ■ ■ = 0, we obtain 


dxleit, 0,0) = ^7 e (*,0,0) = <9*7 e M,0) = 

and so on. Thus u E does not have the C?°°-property near x = 0. 

To prove that the solution is not Q°° across x = 2 1, we consider a point (t, x) near 
this line. We have 


<9 t 7 £ (<,x,0) = —c e (7 e (t,x,0)), a t 2 7e (t,x,0) = c' e (y e (t, x, 0))c E (j E (t, x, 0)), 
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As before, this expression is seen to violate the (/“-property. As e 4- 0; the curve 
x = 7 e (0, 0,t) converges to the line x = 2 1. Hence in every neighborhood of a point 
on x = 2 1, the (/“-property is violated, and so u E is not Q°° across x = 2 1. 

4. Propagation of singularities in wave equations. We study the phe¬ 
nomenon of propagation of singularities in the generalized solution to the Cauchy 
problems for the one-dimensional wave equation with propagation speed depending 
on x 


d 2 u — c{x) 2 d 2 u = 0, t £ [0, T], i£E, 
u| t =o = uo, d t u\ t =o = Ml, 


(4.1) 


and for the wave equation in any space dimension with f-depending propagation speed 


d 2 u — c(t) 2 Au = 0, t £ [0, T], x £ 
u\ t =o = u 0 , d t u\t=o — Mi, x £ 


(4.2) 


In either case, the coefficient c is assumed to be strictly positive and piecewise con¬ 


stant. We work on a finite time interval [0,T), but with arbitrary T > 0. (See the 


comments at the beginning of Subsection 3.1.) 

As outlined in the Introduction, the issue is to find upper and lower bounds on 
the (/“-singular support, that is, to present methods that admit the detection of 
singularities in the generalized solution from its asymptotic behavior in e. 

4.1. The space-dependent case. We begin with problem (4.1). The coefficient 
c(x) is assumed to be given by 


c(x) = c 0 + (ci - cq)H(x) 


where Co, Ci > 0, Cq ^ Ci and H is the Heaviside function. Specifically, we take 


mq = 0, Mi = S(x + 1), 


where <5 is the delta function. We define a generalized function C £ </oo,2(K0 by means 
of a representative (c e ) £S (o,i] = (c * <^ e ) e g(o.i] with a modifier ip e as in Example 1. 
Then singsuppgoo C = {x — 0}. We also produce Ho = 0 and U\ G 1/2, 2 (®0 as the 
class of (<£> e (x + l)) ee (o,i], where ip e again is a modifier as in Example 1. (Actually, 
the modifier need not be the same as the one chosen for the regularization of the 
coefficient c.) Thus we interpret problem (4.1) as the problem 


8 2 U - C 2 d 2 x U = 0 in f/2, 2 ([0, T] x R), 

U\ t =o = U 0 , d t U\ t =o = U 1 ine 2 , 2 (K) 


(4.3) 
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with C, Uq and U\ defined above. By Theorem 3.3, problem (4.3) has a unique solution 
U £ C? 2 , 2 ([0, T] x K.) for any T > 0. We are going to show that, if the coefficient has a 
small jump, then reflection and refraction of the singularity actually do occur at the 
point of discontinuity of the coefficient (see Figure 2). 


t 



Fig. 2 . The Q°°-singular support of the solution U. 

Theorem 4.1. Let cq , c\ > 0 be such that 2 < /co /c\ + \Jc\/co < 4 and let T > 
l/c 0 . Furthermore, let C, Uq and U\ he as described above and let U £ f/2,2([0, T] xR) 
be the solution to problem (4.3). Then it holds that 

singsupp^oo U = {(t,x) | x = — 1 — cgt, 0 < t < T} 

U {(t, x) | x = — 1 + cot, 0 < t < 1/co} 

(4.4) U {(f, x) | x = 1 — cot, 1/co <t <T} 

U {(f, x) | x = -ci/co + c\t, 1 /c 0 <t <T} 

=:r 1 ur 2 ur 3 ur 4 . 


Proof. Assertion (4.4) will hold if we show that 
(4.5) singsuppgoo U C U^ =l r fc . 

To see this, for the sake of simplicity, we assume that T > 2/co- We change the roles 
of t and x and consider {U\ x= -i,d x U\ x= -i) as initial data. By (4.5), the initial data 
(Z7|a; = _i, d x U\ x= ~i) can have singularities only at t = 0 and t = 2/co- We observe 
how their singularities propagate as x goes from —1 to — oo and from —1 to oo. By 
d’Alembert’s formula, U has a representative ((l/2co) J x fi^ot Te(u) ^2/)ee(o,i] near 
(t, x) = (0, —1). Hence singsuppgoo U | x= _i D {0} and d x U\ x =-\ = 0 near t = 0. Let 
us focus on the singularity of 17| a= ^i at t = 0. By d’Alembert’s formula again, it 
propagates along Ti as x goes from —1 to —oo. As x goes from —1 to oo, it propagates 
along T 2 and then splits into a transmitted and a refracted wave at x = 0. This can 
be shown similarly to the proof of Theorem 4.1 of [8]. Thus the singularity of U\ x —-i 
at t = 0 propagates along Ti, T 2 , T 4 and 


r 5 := {(Lz) | x = ci/c 0 - c\t, 0 < t < l/c 0 }. 
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Note that T 5 does not belong to singsuppgoc U. Hence the initial data U | x= _i, 
d x U\ x —-i must have a singularity at t = 2/co which cancels out the singularity on T 5 
emanating from U\ x =~i at t = 0. This singularity at t = 2/co propagates along at 
least T 3 and Ts. It also may propagate along lb. Even if that is the case, it does not 
cancel out the singularity on T 4 emanating from U\ x —~\ at t = 0, i.e., 1?4 belongs to 
singsupp^oo U. In fact, if I ?4 does not belong to singsupp eoo U , then singsupp^oo U 
coincides with the union of Ti, T 2 and IV However, by Proposition 4.3 below, this is 
not possible. Thus (4.5) implies assertion (4.4). 

We next prove that (4.5) holds. Let D be the region bounded by T 3 and V, i.e., 

D := {(t, x) | 1 — cot < x < — C\j cq + ci t, 1/cq < t < T}. 


Indeed, U is zero to the left of Ti and to the right of V U T 4 , and constant in the 
region enclosed by Ti,r 2 ,r 3 . Therefore, 

sing supp^oo U CC (l4 =1 T fe ) U D. 

Thus it remains to show that U is C?°°-regular in D. Put V = dtU — Cd x U and 
W = dtU + Cd x U. Then problem (4.3) can be rewritten as the Cauchy problem for 
a first-order hyperbolic system 


(dt + Cd x )V 

= C'(V - W)/2 

in C? 2 i 2 ([ 0 , T\ x K.) 

(dt - Cd x )W 

= C'(V — W)/2 

in e? 2 , 2 ([0, T] x K.) 

V\ t= o = Vo 

= Ui 

in ^ 2,2 W, 

W\ t = 0 = Wo 

= Vx 

in f/ 2,2 (® 0 - 


We split the initial data (Ui, U\) into two parts, 

{Ui, U\) = (0, U\) + (Ui,Q), 


and let (Vi,Wi) and (IVWV be the solutions to problem (4.6) with initial data 
(0, U\) and (t/i, 0), respectively. Obviously (V), W\) = (0, 0) in D. Hence, if we show 
that (V 2 , W 2 ) = 0 in D , the solution U will be C?°°-regular in D. To do this, we first 
consider the case Co > Ci. The solution (V 2 ,1-V) has a representative (v 2 , e , u ’2,e)ee(o,i]> 
which satisfies the Cauchy problem 


(d t + c e (x)d x )v 2 ,E 

(dt Ce (^X^)d x )'U ? 2,e 

V 2 ,e\t=0 
W 2 , e \t=0 


C' e {x)(v 2 , e - W2, e )/2, 
C' e (x)(v 2 , e - W2,e)/2, 
Ve{%+ 1), 

0, 


0 <t <T, ieK, 
0 <t <T, x G E, 
ieR, 
ieR. 


Consider the characteristic curves 7 + (t, x, r) and (t, x, r) passing through (t, x) at 
time t = t, which are the solutions of the equations 

<9t7 t(t,x,r) = c E (y +(t,x,r)), ^(t,x,t) = x, 

dr%(t,X,T) = -C e (%(t,X,T)), 7 e(t,X,t) = X. 

We can check that 

suppu 2 , e C {(t,x) \ 0<t<T, 

min{—e, 7 +( 0 , -1 -e,t)} <x< 7 +( 0 ,-l + e,t)}, 
suppw 2)£ C {(t,x) \ 0<t<T, 

7 e "(( 1 - 2e)/co,-e,t) <x< min{e, 7 +( 0,-1 + e,t)}}. 
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We now put b e = 1 /c e and (v e ,w e ) = (—b e V 2 t e,b e W 2 } e)y and change the roles of t and 
x. Then (v s ,w E ) satisfies the systems 


(4.7) 


(■ d x + b e (x)d t )v e 
( d x - b e (x)d t )w e 

Ve |x =—e 
w e\x=e 


n e (x)(v e - W E ), 
n e (x){w e -v e ), 
0 , 

0 , 


o <t <t, 

o <t<T, ieR, 
l/co < t < T, 
t e < t < T, 


where 

(4.8) „ M = ±M. = A, logMri 


and t e is such that 7 + (0, —1 + e, t e ) = e. The geometry can be read off from Figure 3. 



Fig. 3. Geometry of the supports of the solution to (4-7)- Note the interchange of the roles of 
x and t. 


The characteristic curves 77 +(i, x, £) and r] e (t,x,£) passing through (t,x) at £ = x 
are the solutions of the equations 


fyr)+(t,x,£) =6 e (0, V t{t,x,x)=t, 

= -b E (0, r)~(t,x,x) =t. 
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Along these characteristic curves, v e and w e are respectively calculated as 
v e (t,x) = -J n £ (£)w s (ri£(t,x,£),£)exp J n e (y)dy^d£, 

6XP (/ ^j 


(4.9) w E (t,x) =w E (r)~(t,x, -e),-e)exp (^J n E (y)dyJ 

-J ne{€)v E (ri~(t,x,£),t)exp J n e (y)dy^dt, 

• exp (/ n e (y)dy 


for 0 < t < T and x £ R such that r]+(t,x, — e) > 1 /cq. Let ii >£ = 7y+(l/co, — e,e) 
and t 2)£ = Vetti ,e; £ ; ~ e )- Then i 2)£ > ti, E > 1/co and i 2>e — 1/co = 0(e) as e 0 (see 
Figure 3). For k £ N, put 


I k = {{t, -e) | l/c 0 + (k- l)(t 2 , E - l/c 0 ) < t < t 2 , e + {k- 1 )(t 2 , E - l/c 0 )} , 
J k = {(£, e) | ti, e + (fc - l)(i 2 , e - l/c 0 ) < t < ti, e + fc(t 2 , £ - 1/co)} ■ 


Define 


(w e ) + (t, -e) = max{0, w E (t, -e)} , 
{w e )-{t,—e) = -min{0,w £ (t, —e)}. 


Let (ui ie ,aci iE ) and (t; 2jE , «; 2i£ ) be the continuous solutions to the integral equations 
(4.9) with initial data (0, (w e ) + (t, —e)) and (0, ( w e )-(t , —e)), respectively. These solu¬ 
tions are obtained by iteration, see [34]. From this, we find that v\ iS (t, x),V2, E (t, x) < 0 
and wi^ e {t,x),W2, E {t,x) > 0 for 0 < t < T and itl such that r]£(t,x,—e) > 1 /cq. 
Using this fact, (4.9) with x = e and (4.8), we obtain, for t > ti )£ , 


(4.10) 


wi, E {t,e) > . — (w E )+(r} e {t,e, -£),-£), 

V ci 


W 2 , E (t,£) > , —(w E )-{r] (t,£,-£),-£). 
V Cl 


Note that w e = wi, e —W2, E and further that w E {t, e) = 0 for t > ti, £ . Hence wi lE (t, e) = 
W2, E {t,£) for t > ti )E . This and (4.10) yield that, for t > ii, £) 


so that 
(4.11) 


W\,s{t,£) = W2, E (t,£) > \j— I w E (r] e (t,£, -e), —e)| 

Cl 


||lCl, £ + W 2 , E 


L“(Ji) 





On the other hand, the solution (t> 3 j£ , w 3tE ) to (4.7) with initial data (0, HiCeHiw/juj,,)) 
at x = — e is given by 




ll w e|U°°(7iUJ 2 )> 


1 

W 3 ,e = g 



[[iCe ||l oo (7 1 U/ 2 ) , 
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and w± tS + W 2 t e is smaller than W 3 e on Ji, i.e., 

( 4 - 12 ) 11^1, e + ^ ( X + l|w e || L =o (JlU72) . 

Combining (4.11) and (4.12), we get 

(4-13) ll^slli 00 ^) - 7 + \f^) ll^e III/ 00 (/1U/2) • 

4 V V c o V c i / 

Since by assumption, 



it must hold from (4.13) that 


||W'e||L“(/iU/ 2 ) = ll^elU“(/i)' 
Together with this, (4.13) leads to 


1 


I Cl 


IKI U-(/ a) <4^- + 

Repeat this process to get for k G N, 


I W e || £■*>(/! )• 


IKIU °°(I k ) < 


k -1 




We now fix to > 1/co arbitrarily and choose k £ G N so that to G /fc e . Then 

1 (to — 1/co)/(^2,£— 1/co) — 1 


IKIUo^j < 





Recall that t 2 , £ — 1/co = O(e) as £ | 0, and note that IliUellioo^) = 0(e p ) for some 
p > 0. Hence, for any N > 0, 

= 0(e N ) as e 4, 0. 

Since t 0 G Ik e and ||'ia e ||x,~(/ fc ) is decreasing in k G N, we get 
||We||L~({(t,-e)|t>t o} ) = 0{e N ) as £ | 0. 

In the region {(t,x) \ t > r]~(to,—e,x), x < —e} the right-hand side of the second 
equation in (4.7) vanishes, so w e is constant along any characteristic curve there, and 
we get 

(4.14) ll“ell£-. ({( t,x ) |t>^- ( t 0 ,- ei * )l a<—e}) = 0 ^ N ) as £ | 0. 

Furthermore, using the fact that v e = 0 on {(£, —e) \ t > to}, we find that, for any 
N> 0, 

ll^llA~({(qx)|t>7 ) +(io,-£,x), x>-e}) = 0( £jV ) as £ 4- 0 , 
U ’z\\L°°({(t,x)\t>ri+(to,-e,x), -s<x<e}) = ) &S £ | 0, 


(4.15) 
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invoking the uniqueness argument as in the proof of [8, Theorem 3.1]. Estimates 
(4.14) and (4.15) are rewritten with ( v e ,w £ ) = {—b e v 2j£ ,b e w 2t£ ) as follows: for any 
N > 0, 


||^2,e|li«.({( tia .)| j 5> toi 7 e -(t 0 ,- £ ,t)<x<-e}) = ^( £ ) aS £ i 0 ) 

11 ^2,e 11 ^oo ({ (t lX ) |t>io, -e<x<'y+(t 0 ,-z,t)}) = as £ | 0, 

ll^ 2 ,ellz,o“({(t,x)|t>t 0 j -e<x<min{ e ,7+(i 0 ,- e ,t)}}) = ^(e N ) aS £ I °. 

This and Lemma 1.2.3 in [20] show that (V 2 ,11' 2 ) = 0 in the region 

{(t, x) \ 1 — c 0 t < x < — ci /cq + cit, t 0 < t < T}. 

Since to > 1/co is arbitrary, we conclude that (V 2 , W 2 ) = 0 in D. 

The case Co < Ci can be treated by the same arguments, using the fact that 
(—r> 2 , £ , W 2 ,s) satisfies the Cauchy problem 


(4.16) 


(dx + b £ (x)dt){-v 2 , e ) 

(d x - b e (x)d t )w 2 ,e 

V2,e\x=-e 
W2,e \x=e 


-n E {x)((-V 2 , e ) + W 2 , E ), 
-n E {x)((-V 2 , e ) + W 2 ,e), 
0, 

0 , 


0 < t < T, id, 
0 < t < T, x G K, 
1/co <t<T, 
t E <t< T, 


and considering the corresponding integral equations for (—v 2tS ,w 2te ). The proof of 
Theorem 4.1 is now complete. 0 

As mentioned in Example 1, one may regularize the piecewise constant propaga¬ 
tion speed c(x) = cq + (c 1 — Co)H(x) in such a way that the corresponding element 
C £ <?oo,2(®0 is l? 00 -regular on R. Indeed, it suffices to use a modifier <ph.( E ), where 
h(e) > 0 and (l//i(e , )) ee (o.i] is a slow scale net. It has been shown in [29] that an 
element C £ £/(R) of bounded type belongs to ^°°(R) if and only if all derivatives are 
of slow scale type. 

Consider problem (4.3) with Uq = 0 and U\ given by the class of (ip e (x+ l)) ee ( 0 .i] 
as in Theorem 4.1. We shall show in the following theorem that the reflected ray 


r 3 = {(t,x) I x = 1 - Cot, 1/co < t < T} 

does not belong to singsupp^oo U , if C is t/°°-regular on R (see Figure 4). In this 
situation, the (^-singularities in the solution arise solely from the flow-out of the 
singularities in the initial data. There are no reflected singularities, just as in the 
classical case of C°°-regular coefficients. This phenomenon has already been observed 
for scalar equations in the ^-setting in [17] and shows again that a non-logarithmic 
regularization is essential for capturing reflected singularities. 

Theorem 4.2. Letco, c± > 0 be such thatco 7 ^ ci and letT > 1/cq. Furthermore, 
let C, Uq and U\ be as described above and let U £ G 2 , 2 ([0,T\ x R) be the solution to 
problem (4.3). If C is Q°°-regular on R, then 

singsuppgoo U = {(t, x) | x = —1 — cot , 0 < t < T} 

U {(t, a;) | x = — 1 + cot, 0 < t < 1/co} 

U {(t, x) | x = — ci/co + c\t, l/c 0 <t<T} 

=-. Ti ur 2 ur 4 . 


(4.17) 
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t 



Fig. 4. The Q°°-singular support of the solution U with slow scale coefficient. 


Proof. Assertion (4.17) will hold if we show that 
(4.18) singsupp go o [/C Ti U r 2 U T 4 . 

This can be seen as follows. For the sake of simplicity, we assume that T > 2/cq. 
As in the proof of Theorem 4.1, we interchange the roles of t and x and consider 
[U\ x =-\,d x U\ x =-\) as initial data. From (4.18) and d’Alembert’s formula, we see 
that singsupp^oo U\ x =-i = {0} and d x U\ x =-i = 0 near t = 0. The singularity of 
U\ x --\ at t = 0 propagates along Ti as x goes from —1 to — oo. As x goes from —1 
to oo, it propagates along r 2 and No splitting of the singularity occurs at x = 0 
unlike in Theorem 4.1. This can be proven similarly to the proofs of Theorems 4.1 
and 5.1 of [8]. Thus (4.18) implies assertion (4.17). 

To show that (4.18) holds, we consider system (4.6). Put Fg = {(t, 0) | 1/co < 
t < T}. Let D\ be the region bounded by F 3 and Tg, and U 2 be the region bounded 
by T 4 and Tg. It is immediate to check that 


suppP c r 2 u r 4 u r 6 u n 2 , 
supp w c Ti u r 3 u r 6 u D\. 


Put B = 1/C and (V, W) = (—BV,BW). Clearly 


(4.19) 


suppP c r 2 u r 4 u r 6 u d 2 , 
supp w c Ti u r 3 u r 6 u D\. 


Furthermore, (V, W) satisfies the system 

(d x + Bdt)V = N{V-W) % 
(8 X - Bdt)W = N(W - V), 


where N = B'/(2B). Take the value of (P, W) along some line {x = zo} with xq > 0 
as initial data and consider the system above. Using the commutator argument as 
in the proof of Theorem 5.1 of [8], the derivatives ( d x + Bdt) k W can be estimated 
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in terms of the H°°-norm of the zeroth derivative of V for whatever power k. The 
derivatives ( d x — Bd t ) k W are estimated on Di using the fact that V vanishes there. 
We conclude that W is (^““-regular on 1?3 U D\. With this in mind, consider now the 
value of (H, W) along {x = Xi} with X\ < 0 as initial data. Then Proposition 4.6 
below shows that W is C?°°-regular along To and V is t/°°-regular on To U £> 2 - Hence 
from (4.19), 


singsuppgoo HcT 2 UT 4 , 
singsupp^oo W C Ti. 

This implies (4.18). The proof of Theorem 4.2 is now complete. 0 

4.2. The time-dependent case in one dimension. This subsection is de¬ 
voted to problem (4.2) in one dimension: 


(4.20) 


d 2 u — c{t) 2 d 2 u = 0, t > 0, x £ R, 
u|t =0 = y-o, d t u\t =o = Ui, X £ R. 


The case of delta functions as initial data has been settled in [8]. The goal of this 
subsection is to extend the results of [8] to arbitrary C?°°-point singularities in the 
initial data. Assume that the coefficient c(t) is given by 

c(f) = c 0 + (ci - c 0 )H(t - 1) 


where Co, c 4 > 0, Co ^ ci and H is the Heaviside function. As in the previous 
subsection, we define a generalized function C £ (7(R) by means of a representative 
( c e)ee(o,il = (c* ¥>e) e g(o,i] with a mollifier ip e as in Example 1. Then singsupp^oo C = 
{f = 1}. Thus we will consider the Cauchy problem 

d?U -C 2 d 2 U = 0 in£?([0,oo)xR), 

U\t=o = t/ 0 , dtU | t=0 = Hi in g(R). 

By [8, Theorem 3.1], problem (4.21) has a unique solution H £ 5([0,oo) x R) for any 
initial data Ho, Hi £ £(R). 

An important ingredient in the arguments is the persistence of regularity from 
regular initial data. Thus we are first going to show that, if Ho, Hi £ £7°°(R), then H 
is f?°°-regular off the t/°°-singular support of the coefficient C. 

Proposition 4.3. Let C be as described above and let U £ £/([0,oo) x R) be the 
solution to problem (4.21) with initial data Ho, Hi £ C?(R). If Ho, Hi £ t/°°(R), then 

(4.22) singsuppgoo H C {(f, x) \ t = 1, x £ R} . 


Proof. The proof is similar to that of [8, Theorem 3.1], where an existence result 
for problem (4.21) with a more general coefficient has been proven. Put V = dtU — 
Cd x U and W = dtU + Cd x U. Then problem (4.21) can be rewritten as the Cauchy 


problem for 

a first-order hyperbolic system 




(d t + Cd x )V 

= M(V - W) 

in 0([O, oo) 

x R) 

(4.23) 

(d t - Cd x )W 

= M{W — V) 

in S([0, °°) 

x R) 

V\t=o = Ho 

= Hi - (C7| t= o)H5 

in 1/(R), 



VH| t =o = W 0 

= Hi + (C|i = o)Ho 

in fy(R), 
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where M = C'/(2C) £ Cy(R). The solution (V,W) has a representative (v e , w E ) se ( 0tl ] 
satisfying the Cauchy problem 


(d t + c E (t)d x )v E = fi E (t)(v E -w E ), t> 0 , iel, 

(d t - c E (t)d x )w E = /j, E (t){w E -v E ), t>0, 

Ve\t=0 = v 0e = Ule - C e (0)u' Oe , iel, 

W E \ t =0 = WOe = Ui E + C E (0)u' Oe , l£l, 


(4.24) 


where (uo e ) ee (o,i]> ( M ie)ee(o,i], (c E ) Ee ( 0 ,i] and (Me)ee(o.i] are representatives of U 0 , U\, 
C and M, respectively, such that n e = d E /{2c E ). Consider the characteristic curves 
7 +(i, x , r) and 7 ~(t, x, r) passing through (t, x) at time r = t which satisfy 


d T r= c e (t), 7 £(t,x,t)=x, 

d T % (t,x,r) = -c e (t), % (t, x, t) = x. 


Along these characteristic curves, v E and w E are respectively calculated as 



(4.25) 


For each T > 0, we define Kt as the trapezoidal region with corners (0, —£), (T, + 

C 2 T), (T, £ — C 2 T), (0,^), where C 2 = max(co,ci). From (4.25), the following inequal¬ 
ities are easily deduced: 



We add these two inequalities and apply Gronwall’s inequality to get 


v e\\l°°(K t ) + ll^ell L°°(K t ) 



Differentiating the equations and using the same argument, we obtain for any a £ No, 


ll<9> E || l°°{k t ) + \\d2w E \\ L ao( KT) 



(4.26) 


On the right-hand side, the terms involving vq e and wq s are of order 0(e N ) for 


some N > 0 independent of a. The exponential term is uniformly bounded in e. 
Note that, on every compact set outside {(t,x) \ t = 1, x £ R}, d t v e = —cd x v E and 
dtW e = cd x w e , when e is small enough. This and inequality (4.26) imply that V and 
W are £/°°-regular outside \ t = 1, x £ R}. From the definitions of V and W, 

assertion (4.22) follows. 0 
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The central issue of this subsection is the propagation of the C?°°-singular support 
from arbitrary initial data in £7(R). As the prototypical case, we focus on a point 
singularity, say at the origin. We take Uq, U\ £ £°°(K \ {0}). Choose a cut-off 
function \ £ Co° (®0 which is identically equal to one in a neighborhood of the origin. 
Split the initial data (Uq, U\) into two parts, 


(f/o, U 1 ) = ( X Uo, X Ui) + ((1 - X)U 0 , (1 - x)Ui). 


Then the t/^-singular support of the solution U to problem (4.21) with initial data 
(yt/o, XUi) is contained in the set 


(4.27) 


•|(t, x) | x = ± J c(s) ds , t > o| 

U |(t, x) | x = ± ^2 J c(s) ds — J c(s) ds^j , t > 11 


in the limit as the support of x shrinks to {0}. This is most easily seen by introducing 
V = d t U — Cd x U, W = d t U + Cd x U, which satisfy a system of equations similar to 
(4.23) with a corresponding cut-off in the initial data and hence are supported in a 
neighborhood of the lines defined by (4.27). Furthermore, by Proposition 4.3, the 
solution to problem (4.21) with initial data ((1 — x)Uq, (1 — X)U\) is £°°-regular off 
{(t,x) | t = 1, x & R} for any x as above, since (1 — x)Uq and (1 — xW'\ belong to 
t7°°(R) for any x as above. We summarize the result in the following proposition (see 
also Figure 5). 

Proposition 4.4. Let C, Uo and Ui be as described above and let U £ f?([0, oo) x 
R) be the solution to problem (4.21). Then the set 


(4.28) 


{(t, x) | t = 1, x £ M} U | (t, x) | x = ± J c(s) ds, t > 0 j 
U (t, x) | x = ± ^2 J c(s) ds — J c(s) ds ^ , t > 11 


is an upper bound for the Q°° -singular support of the solution U to problem (4.21). 
In the case of delta function initial data, we have shown in [8] that the t7°°-singular 
support of the generalized solution coincides with the set (4.27). The line {t = 1} 
generally may belong to the t7°°-singular support of the solution, even for C°°-initial 
data, as shown by the example 

. . x 2 

u e {t,x ) = — + 

which solves the wave equation dfu e — c e (t) 2 df.u e = 0 with initial data u £ ( 0, x) = x 2 /2, 
dtu s (0,x) = 0. In general, the t7°°-singular support of the generalized solution may 
be a proper subset of the set given by (4.27). This depends on the interplay between 
the initial data and the coefficient c(t). 

In the following theorem, we give a necessary and sufficient condition under which 
initial singularities propagate out and the transmitted rays 



r+ : 
r_ : 


(t, x) | x = J c(s) ds, t > 01 , 
|(f, x) | x = — J c(s) ds, t > 01 
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belong to singsuppgoo U. 


t 



Fig. 5. An upper bound for the Q°° -singular support of the solution U. 


Theorem 4.5. Let C, Uq and U\ be as described above and let U E £/([0, oo) x R) 
be the solution to problem (4.21). Then: 

(a) */ sing suppgoo (U\ - (C\ t =o)U' 0 ) = {0}, then 

r + = | (t, x) | x = J c(s) ds, t > o) C singsuppgoo f/; 


(b) */ sing suppgoo (J7i — (C'|t=o)17o) = 0, thenU is Q°°-regular along r+\{(l,co)}; 

(c) if sing supp g =o (t/i + (C\ t =o)Uo) = {0}, then 

r_ = j(t,x) x = — J c(s) ds, t > o| C sing suppgoo U; 


(d) */sing suppg«, (i7i+(C'| t= o)C/o) = 0, thenU is Q°° -regular along T_\{(1, — Co)}. 
Proof To show that (a) and (b) hold, we consider the solution (V, W) to problem 
(4.23) with (Vo, Wo) = (x(f7i — (C|t=o)I7o), 0)- As may be seen from the proof of 
Proposition 4.3, if sing suppgoo {U\ — (C'|t = o)E/’o) = then ^ is ^“-regular outside 
{{t,x) | t = 1, i£l} for any y. From this, case (b) follows. As for (a), it suffices to 
show that if sing suppgoo {U\ — (C'|t=o)f7o) = {0}, then 


(4.29) 


sing suppgoo V = T + 


holds in the limit as the support of x shrinks to {0}. It is easy to see that supp kcT+ 
holds as the support of x shrinks to {0}. Hence, if we show that for any x, V is not Q°°- 
regular on R for any t > 0, then (4.29) will follow. We take (v E , ty e )e£(o,i] satisfying 
problem (4.24) with (uo e , wq e ) = (x(ui e ~ c e(0)u' Oe ), 0). Differentiate the first equation 
of system (4.24) in x, k times, multiply by 2(d 1 fv e ) and integrate to obtain 


d_ 

dt 


(d£v E ) 2 dx 


2 He(t) / (d^v e ) 2 dx - 2fj, e (t) / (d%v £ )(d%w e )dx. 

J—oo J—oo 


Similarly 


d 

dt 


(d%w E ) 2 dx 


2ti E {t) / {d*w e ) 2 dx - 2n E (t) / (d%v E )(d%w £ )dx. 
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Take their difference and integrate to get 

/ OO 

{(dxVe) 2 ~ (d*w £ ) 2 ) dx 

-CO 

= J ({dxVoe ) 2 - (d^w 0e ) 2 )dx ■ exp (2 J /j, £ (s) ds'j . 

Since Wo e = 0 and = c' e /{2c e ), it follows that 

(4.30) f (d*v e ) 2 dx > f (d*v 0e ) 2 dx. 

J —co ^e(O) J _oo 

Note that vq £ vanishes outside some compact set, independently of e G (0,1]. Then 
it is clear from finite propagation speed that, for some constant C\ > 0, 

/ CO 

(d k x v e ) 2 dx < c^dSWl-m- 

-CO 

The Sobolev embedding theorem yields that, for some constant Ci > 0, 

(4-32) IldSellW) < C2 E H^ll t 2 (R)' 

(<k +1 

Combining (4.30)-(4.32), we get, for some constant C 3 > 0, 

(4-33) ll^x ?; 0e|||oo(R) < C 3 ll^f t ’e|li°°(R). 

e<k +1 

By assumption, Vo = x{U\ — {C\ t =o)U'o) is not 0°°-regular on R for any %. This and 
inequality (4.33) imply that for any y, V is not (/“-regular on R for any t > 0. Thus 
(4.29) follows. Cases (c) and (d) can be argued similarly. □ 

When the coefficient C in problem (4.21) is (/“-regular, the same phenomenon 
as in Subsection 4.1 occurs: singularities can only be produced by the initial data, 
and there are no refracted rays, as will be seen next. Again, we use a mollifier ^Ph(e)i 
where h(e) > 0 and (l/A(e)) e e(o.i] a s l° w scale net to regularize c(t). Then the 
corresponding element C £ (/(R) belongs to (/“(R). First, we have the following 
global regularity result, which can be proven similarly to the proof of Proposition 4.3. 

Proposition 4.6. Let C be as described above and let U £ <7([0, 00 ) x R) be the 
solution to problem (4.21) with initial data Uq, U\ £ (/(R). If C, Uq, U\ £ (/“(R), 
then U £ £?“([0,oo) x R). 

Turning to propagation of a point singularity, we take Uq, U\ £ (/“(R \ {0}) as 
in Theorem 4.5. Again, the cut-off argument employed in the derivation of Propo¬ 
sition 4.4 yields that the (/“-singular support of the solution U to problem (4.21) is 
contained in the set given by (4.27), noting that no singularity is present along the 
line {t = 1}. 

Next, the same argument as in the proof of Proposition 5.1 of [8] shows that the 
refracted rays 


(t, x) | x = ± (2 J c(s ) ds — J c(s) ds ^ , t > l]. 
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do not belong to the C/°°-singular support of the solution U , and so the singularities 
can occur only along the transmitted rays (see Figure 6) 

(4.34) | (t, x) | x = ± J c(s) ds , t > 0 

Proposition 4.7. Let C, Uq and U\ be as described above and let U £ C?([0, oo) x 
R) be the solution to problem (4.21). Then the set (4.34) is an upper bound for the Q°°- 
singular support of the solution U to problem (4.21). As for the question whether the 
singularities actually propagate out, the same argument as in the proof of Theorem 
4.5 can be applied and thus the following theorem holds. 


t 



Fig. 6. An upper bound for the Q°°-singular support of the solution U with slow scale coefficient. 

Theorem 4.8. Let C, Uo and U\ be as described above and let U € !?([0, oo) x R) 
be the solution to problem (4.21). Let C £ £7°°(R). Then: 

(a) if sing suppgoo (U \ - (Cl^o)^) = {0}, then 

T + = |(f, x) | x = J c(s) ds, t > 0 j> C singsuppgoo [/; 

(b) z/singsupp(Lb — ( C\t=o)U'o ) = 0, then U is Q°°-regular along T + ; 

(c) if sing supp a .o (C/i + {C\t=o)U'o) = {0}, then 

T_ = | (f, x) | x = — J c(s) ds, t > o| C sing supp^oo U ; 

(d) z/singsuppgoo(C/i + (C\ t =o)Ub) = 0, then U is Q°°-regular along T_. 

4.3. The time-dependent case in multiple dimensions. We here consider 
problem (4.2) in multiple dimensions: 

dfu — c(t) 2 Au = 0, t £ [0, T], x £ R d , 

u|t=o = uo, d t u\ t =o =u\, x £ R d . 

As in the one-dimensional case, 

c{t) = c 0 + (ci - co)H(t - 1), 

where cq , C\ > 0, cq ^ C\ and H is the Heaviside function. Specifically we take 


u 0 = 0, 


Ml = S, 
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where <5 is the delta function. We define a generalized function C £ 0oo,2[O, oo) by 
means of a representative (c e ) ee ( 0 ,i] = ( c * V 3 £)ee(o,i] with a modifier p e as in Example 
1. We also produce Uo = 0 and U\ £ t/ 2 , 2 (R d ) as the class of (v? e (|®|))ee(o,i], where 
ip e again is a modifier as in Example 1. (Actually, the modifier need not be the same 
as the one chosen for the regularization of the coefficient c.) Thus we will consider 

d 2 U - C 2 AU = 0 in C?2, 2 ([0, T] x R d ), 

(4.35) 

U\t=o = U o = 0, dtU\ t =o = U 1 in0 2 , 2 (R d ). 


By Theorem 3.5, problem (4.35) has a unique solution U £ (7 2 , 2 ([0,T] x U d ) for any 
T > 0. 

We begin with the odd-dimensional case. It was shown in [8, Theorem 4.1] that, 
in the one-dimensional case (d = 1), it holds that 


(4.36) 


singsupp^oo U = |(£, x) | |x| = J c(s)ds, 0 < t < t| 

U / (t, x) | |x| = 2 J c(s)ds — J c(s)ds 


, 1 <t <T > . 


We are going to show that the same holds for ad odd d > 3. 

Theorem 4.9. Let d > 3 be odd. Assume that C, Uq and U\ are as described 
above and further that U £ £/2,2([0,T] x R d ) is the solution to problem (4.35). Then 
the Q°°-singular support of U coincides with the set given by (4.36). The support of 
U also coincides with the set given by (4.36). 

Proof On the level of representatives, problem (4.35) reads 


dfu e — c e {t) 2 Au E =0, 0 < t < T, x £ R d , 

«e|t=O = 0, d t U e \t=o = Ve{\x\), X £ R d . 


Since the initial data are radially symmetric, so is the solution u E . Put Ud, e {t,r) := 
u e (t,x) for |r| = |x|, where the subscript d represents the dependence of Ud, E on 
dimension. Then Ud,e satisfies the Cauchy problem 

( 4 37 j dt u d,e ~ C e {t ) 2 (d 2 + -<9 r ^ u di e = 0, 0 < t < T, r £ R, 

'U J d,e\t=0 b, fdt'U’d,e\t=Q Te: ^ 

According to [19], if Ud, e satisfies the first equation in (4.37), then (—1 /r)d r Ud, E sat¬ 
isfies the first equation of (4.37) with d replaced by d + 2. This yields that, if v n , E , 
n £ N, satisfies the Cauchy problem 

d?v n , E - c e (t) 2 dfv n , e =0, 0 < t < T, r £ R, 

^n,e|t=0 — 0) V £ R, 

d t V n ,e\t =0 = flooi-Sn) ' ' ' f^ 00 (-S 1 )(p E (s 1 ) ds i • • • ds n , T £ R, 

then it 2n +i,e : = [(—1 /r)d r ] n v ntE satisfies problem (4.37) with d = 2n + 1. It is clear 
from the definition of <p E that dtv n:E \t=o is symmetric, (dtV ntE \t=o)' > 0 on [—e, 0] 
and supp<9tU n , £ |t=o C [— e,e\. Hence the same argument as in the proof of Theorem 
4.1 of [8] shows that the ^“-singular support of the class of (i> n ,e)ee(o,i] and the 
support of the class of (9 r u n ,e)ee(o,i] are equal to the set given by (4.36). This and 
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Fig. 7. An upper bound for the Q°°-singular support of the solution U. 


the relationship U 2 n +i,e = [(—1/ r)d r ] n v n:E imply that the C?°°-singular support and 
the support of the class of (w 2 n+i,e)ee(o.i] are the same as the set given by (4.36). 
Thus the assertions follow. □ 

In the even-dimensional case, the following result holds (see Figure 7). 

Theorem 4.10. Let d be even. Assume that C, Uq and U\ are as in Theorem 
4.9 and further that U £ f?2,2([0,T] x K d ) is the solution to problem (4.35). Then it 
holds that 


(4.38) 


singsupp^oo U C |(f, x) | t = 1, |x| < J c(s)ds| 

U | (t, x) | |a;| = J c(s) ds , 0 < f < T j- 
U / (t, x) | \x\ = 2 f c(s)ds — J c(s) ds 


, 1 <t<T 


Proof. As in the proof of Theorem 4.9, we consider problem (4.37). If v ntE , n £ N, 
satisfies the Cauchy problem 

dt v n,e - c e {t ) 2 (df + v n>s = 0, 0 < t < T, r £ R, 

'Vn,e\t=0 — 0? T G K., 

(4.39) d t v n ^\ t=0 

{ f^ 00 (-s n -i)---fl 2 00 (-Si)ip e (s 1 )ds 1 ---ds n -i, r £ R, if n > 2, 

\ Te(r), r £ R, if n = 1, 

then U 2 n,e ■= [(—1 /r)d r ] n ~ 1 v niS satisfies problem (4.37) with d = 2 n. By the argu¬ 
ment as in [39, p.66], the solution v n>s of problem (4.39) is given by 


(4.40) 


f 1 1 

v n ,e{t,r)= w n , £ (t,rp)dp, 

Jo y/1 - p 2 
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provided that (ui n ,e)e€(o.i] £ £m([ 0,T] x R) is symmetric and satisfies 

(4.41) dtW ntS - c e (t) 2 d 2 w nt£ = 0, 0 < t < T, r e R, 

^n,e\t—0 0 , I £ R 5 

f 1 1 

(4.42) d t v„ }e \t=o= / , (w„ tS ) t {0,rp)dp, r G R. 

Jo 

(For clarity of notation, we write (w n:E )t for (dtw n ^ £ ).) The existence of such w n>£ can 
be seen as follows. By the change of variable y = rp , we obtain from (4.42) that 


r Irl /' |r| 1 

d t v n , E \ t=0 = — . (w n , E ) t {0,y)dy= > „, E ) t (0, y) dy. 

Jo r^/r z -y z Jo \/r z -y z 


V r2 ~ 


Thanks to [38], this integral equation can be solved and <9tU> n , e |t=o is given by 


(4.43) 


0 — 


ld_ 

it dr 


2|y| p 

V 1 - p‘ 


■■{v n , e )t{.0,rp) dp, 


so that (9tru„ )e |t-o) ee (o l i] belongs to Em (R). This guarantees that (w) n ,e)eg(o,i] exists. 
Note that the integrand is a smooth function of r, due to the symmetry of (u n , e )t(0, •). 

We now investigate the regularity of (itfn,e)ee(o,ii- Since supp dtV UtE \t=o C [— e, e], 
we see from (4.43) that (dtWn,e|t=o)ee(o,i] belongs to £^f(R \ {0}). Therefore, as in 
the previous subsection, we can show that (w„ i£ ) £g (o,i] belongs to £|f([0,T] xR\T) 
with T given by 


T := {(t, r) 11 = 1, r£ R} U |(f,r) | |r| = J c(s ) ds, 0 < t < T 


U (t,r) | |r| = 


2 / c(s) ds c(s) ds 
Jo Jo 


1 <t<T\. 


We next show that (un, £ ) eg (o.il also belongs to £^J([0, T] x R\T). To this end, we 
take a cut-off function x £ Co° (®0 which is identically equal to one in a neighborhood 
of the origin. Split dtw nt£ \t=o into two parts 


dtWn,e\t=o = x(r)d t w n , E \t= 0 + (1 - x( r ))dtW n ,s\t=o- 


For anyx, (x3tWn, e |t=o) ee (o,i] £ ^( R \{°}) and (( 1 -x)9t« , n, E |t=o) £ 6(o,i] £ £m( r )- 
Let w n ^ i iE and w ra , 2 ,E be the solutions corresponding to the data x9tW n , E \t=o and 
(1 — x)dtWn,e\t=o, respectively. From the linearity of equation (4.41), w n<E = w nt i iE + 
w n ,2,E- By (4.40) we have 

v n ,e(t,r) = [ 1 w„,i tE (t,rp)dp+ f w n ^, E {t,rp) dp 

Jo yj\-p z Jo \J\ — P z 

— . I n ,l,e{t,r) + I n , 2 ,e(t,r). 

It is clear from Proposition 4.3 that (wn,2,e) £g (o,i] £ ^M({(*! r ) I t ± 1, r £ R}) 
for any x- This implies that (/„,2, E )ee(o,i] £ ££?({(*> r ) I t ± 1, r £ R}) for any x- 
Hence, in order to show that (u„ )£ ) eg ( 0i i] £ £^([0, T] x R\T), it suffices to prove that 
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(4i,l,e)ee(0,l] S £}% ([0, T] x R \ r) in the limit as the support of x shrinks to {0}. We 
divide [0, T\ xl\T into the following six regions (see Figure 8 ): 

D\ := |(t,r) | |r| < J c(s)ds, 0 < t < l| , 

D 2 := | |r| < 2 J c(s) ds — J c(s) ds, 1 < t < 1 + co/ci| , 

D 3 := |(i,r) | |r| < J c(s ) ds - 2 J c(s) ds, 1 + co/ci < t < t| , 


Da := < (t, r) 


2 f c(s) ds — f c(s) ds < |r| < f c(s) ds, 1 < t < t! , 

Jo Jo Jo J 


D 5 := j(t,r) | |r| > J c(s) ds, 0 < t < l|, 
D 6 ■.= | (t, r) | |r| > J c(s) ds, 1 < t < t| . 


It is easily checked that (4i,i,e)ee(o,i] £ Dj) for any y, using the fact that 


t 



Fig. 8. The singular support in even space dimensions, radial cross section. 

(w ; n,i,e) E e(o,i] G £m(Fj =1 Dj) for any %. As for the other regions, we first consider 
Da D {(t,r) | t < 1 + Co/c\, r > 0}. Let K C Da fl {(t,r) | t < 1 + cq/ci, r > 0} 
be any fixed compact set. We start by proving the following support property of the 
derivatives of w n .\. e : 

(SP) For any (t,r) £ K, there exist 0 < r 0 , tt < r independently of e, such that 
(w n ,i,e)r{t,r') = (wn,i, E ) t (t, r') = 0 for 0 < r' < To and n < r' < r, possibly 
after shrinking the support of %■ 

Indeed, w Ut i )£ solves the problem 

d/w n ,l,e - C e {t) 2 d/w n ,i,e = 

' l ^n.l : e\t=0 = 

StW n ,l,e\t=0 = 


0 , 

0 , 

XdtW n ,e\t= 0 , 


t > 0, r £ R, 
r £ R, 
r £ R, 
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SO 


fn,£ — ifJt l,e? 

gn, e = (dt + C e dr)w n ,\,e 


satisfy the corresponding first order hyperbolic system 


(d t + C e d r )f nt£ = - 9 n,e), t > 0, r £ R, 

(^t ^£^r)Sn,£ = /n,e)? t > 0, I E R, 

/n,e|t=0 = X'9tWn, E |t=0, r G R, 

5n,e|t=0 = XdtW n ,£\t=0, rGR, 

where /z e (t) = c' e (t)/2c e (t). The supports of f n>e and g n , e are contained in the set 
•| (f , r) | r = ± J c(s) ds , t > 0 1 

U (t, r) | r = ± ^2 J c(s) ds — J c(s) ds^j , t > 11 


in the limit as the support of x shrinks to {0}. It follows from the definitions of /„ j£ 
and g n , E that the supports of d t w n i_ £ and d r w n p, E are contained in the same set as 
the support of x shrinks to {0}. This proves the indicated support property (SP). 

We further observe that c E (t) = Ci on K for s £ (0,1] small enough. Taking 
the 2m-th derivative of 1>e with respect to t and integrating by parts, we get, for 
(t, r) £ K and e £ (0,1] small enough, 


2m /*! 1 

df m I nXe {t,r) = / df n - 1 (w nXE ) r (t,rp)dp 

r Jo v 1 — P 


(- 1 ) 


2 m— 1 „ 2 m 




{WnXe)r{t,Tp) dp. 


Note that the integration by parts is justified, because (w n p :E ) r (t,rp) = 0 near p = 0 
and p = 1, due to property (SP). Similarly 


d? m ~ 1 I nA , E (t,r) = 


(- 1 ) 


2m—2 Jm-2 



2m—2 


v/W 5 


K,i ,£)t(t,rp) dp. 


Hence, for any a £ No, there exists k > 0 such that 


||<9 t “/„,i, e ||£°o(ir) < ftmax{||u) n! i i£ || i oo([ 0)T ] xR ), ||3 r uvi,e||L°°([o,T]xR)> 

||0tW n ,l,s||.L°°([O,T]xR)} 

for sufficiently small e £ (0,1]. Analogous estimates for the r- and mixed derivatives 
can be obtained. The same argument applies in the case that K C (D4 U D 5 U D 6 ) f~l 
{(£, r) | f ^ 1 + co/ci, r £ R} is any fixed compact set. Therefore, taking into account 
finite propagation speed, we obtain that (/ n ,i,e)ee(o,i] 6 4 U D 5 U D 6 ) as the 

support of x shrinks to {0}. Thus (I n , i,e)ee(o,i] £ f^([0,T] x R\T) as the support 
of x shrinks to {0}. 
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We finally prove that (w2n,e)ee(o,i] €= £^([0, 7 1 ] x R \ T') with T' given by 

r' := j(i,r) \t = l, |r| < J c(s)dsj 

U | (f, r) | |r| = J c(s ) ds, 0 < t < T | 


U (t,r) | \r\ = 


2 / c(s) ds c(s) ds 
Jo Jo 


1 < t < T 


This means that assertion (4.38) holds. From the first equation of problem (4.39), we 
have 

s d?v n , e (t,r) -c £ (f) 2 <9 2 u n , £ (t,r) 

— d r v n , E (t,r ) =-—2--. 

r c e {ty 


Note that, on every compact subset outside {(f, r) \ t = 1, r £ R}, c e {t) = c, when 
e > 0 is small enough. Therefore ((—l/r)9 r .u rai£ ) ee (o i i] £ £^([0,T] x R\T). A direct 
calculation shows that (—1 /r)d r v n , e satisfies the first equation in problem (4.37) with 
d = 4. Thus we get 



<9 t 2 ((-l /r)d r v n ^ - c e (t) 2 d^({~l/r)d r v n} ^j 
3 c e {t) 2 


so ([(-l/r)5 r ] 2 u ni£ ) eg(0jl ] £ f^([0,T]xR\T). By induction, we see that (w 2 n, £ ) eG (o,i] = 
([(—l/r)9 r ] n_1 u n)£ ) £g ( 0) i] £ £^([0,T] x R \ T) for any n £ N. In addition, by finite 
propagation speed, 


( u 2 


: n , e ) ee( o,i]£^({(t,r)||r|> J c(s) ds, 0 <t <T . 


Thus (u2n,e) eG (o,i] £ x R \ T') for any n £ N. The proof of Theorem 4.10 is 

now complete. □ 

Remark 8. If the coefficient c(t) = Co + (ci — co)H(t — 1) is regularized in such a 
way that the corresponding element C £ f7oo,2(R) is f?°°-regular on R, then as in the 
one-dimensional case, one can show that 


singsuppgc*, U 


(t, x) | |ir| = J c(s ) ds, 0 < t < 


in any space dimension. 

5. Appendix: Associated distributions. In case the initial data and the 
coefficients are sufficiently smooth, hyperbolic equations and systems possess classical 
or distributional solutions. In numerous such situations it has been verified that the 
Colombeau solution admits the classical solution as associated distribution. This is 
generally the case when the coefficients are C°°-smooth (see e.g. [31, 34]). An instance 
with merely L°°-coefficients is the linearized Euler system of isentropic gas dynamics, 
in which the association result has been proven in [18, 33]. For the one dimensional 
wave equation with a coefficient depending on time only and suffering a jump, i.e., 
problem (4.20), the Colombeau solution was shown to be associated with the piecewise 
distributional solution for arbitrary distributions as initial data in [8]. 
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If the coefficients are piecewise constant with a jump across a smooth hyperplane 
in space, the problem may be considered classically as a transmission problem. The 
one-dimensional problem 

(5.1) a(x)d 2 z(t,x) — d x (b(x)d x z(t, x)) = 0 
with initial conditions 

(5.2) z( 0, x) = z 0 (x), d t z( 0, x) = zx(x) 

has been treated in [18], assuming that a and b are strictly positive and piecewise 
constant with a jump at x = 0. Viewing (5.1) as a transmission problem across 
x = 0, and assuming that zo is continuously differentiable and Z\ is continuous, 
there is a unique distributional solution ~z £ 2V(R 2 ), the classical connected solution , 
characterized by the following properties: 

- z is a distributional solution to (5.1) on the open half planes {x < 0} and 
{x > 0}; 

- z £ C([0,oo) x R) and z( 0, •) = zq; 

- -j^z £ C([0, oo) x R) and -§fz(0, •) = z i; 

- the function x —>• b(x)-^z(t,x) is continuous for almost all t> 0. 

A Colombeau solution Z £ Q (R 2 ) can be constructed corresponding to initial data 
i(zo), a(2i), with suitable interpretations of a(x), b(x) as elements A, B of f/(R). 
It was shown in [18] that the Colombeau solution Z £ fy(R 2 ) admits the classical 
connected solution ~z as associated distribution. The problem (3.9) is a special case of 
(5.1). 

It is the purpose of this appendix to compute the classical connected solution, to 
extend the calculation to delta functions as initial data, to show that the Colombeau 
solution is associated with the distributional solution in this case as well, and to verify 
that the location of the singular support of the associated distribution coincides with 
the C7°°-singular support computed in Theorem 4.1. We begin by computing the 
classical connected solution to the transmission problem 

d 2 u(t 7 x) — c(x) 2 d 2 u(t , x) = 0, t > 0, x £ R, x ^ 0, 

(5.3) 

u(0,x) = uq(x), dtu(0, x) = ui(x), iff 

where 


c(x) 


c_, x < 0, 
c+, x > 0 


and both c_ and c+ are strictly positive. As noted above, the classical transmission 
conditions at x = 0 are 

- u is continuous across x = 0 and 

- d x u is continuous across x = 0. 

The first condition entails that also dtu is continuous across x = 0. We introduce 
v = (d t - cd x )u : w = (d t + cd x )u : 


vq(x) = ui(x) - c(x)u' 0 (x), 


wq(x) = u\{x) + c(x)u' 0 (x) 
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and rewrite equation (5.3) on either side of the line x = 0 as the system 

(dt + cd x )v = 0 
(5.4) (d t - cd x )w = 0 

11 ( 0 , 2 ) = 20 ( 2 ), 12 ( 0 , 2 ) = 700 ( 2 ) 


where c = c+ for 2 > 0 and c = c_ for 2 < 0. Note that we can recover 


d t u = -(v + to), 


d x u = —(to 
2 c v 


v), 


u(t, 2 ) = 2 - 0 ( 2 ) + / dtu(s,x)ds. 

Jo 

We shall use the suffix ± to indicate parts of the solution on 2 < 0 and on 2 > 0. 
Consider the regions displayed in Figure 9, which are separated by the lines 2 = —C—t, 
x = 0, 2 = c+t (plotted with c_ = 1, c+ = 2). 


t 



Fig. 9. Regions of the classical transmission problem. 


Clearly, 

v-(t,x) = vo(x — c-t) in regions (I) and (II), 
v+(t,x) = 00(2 — c+i) in region (IV), 

70 +(i, 2 ) = 700(2 + c+t) in regions (III) and (IV), 

70_(f, 2 ) = 100(2 + C-t) in region (I). 

It remains to determine 0 + in region (III) and io_ in region (II). To this end we recall 
that |(o + 10 ) and ^(10 — 0 ) should not jump across 2 = 0. This yields 


0 + + 10+ = V- + W- 
^( 10 + - 0+) = (lO— - 0_) 

along 2 = 0. Solving this linear system of equations results in 

2c+ c_ — c+ 2c_ c+ — c_ 

0 + =- V- H- 10 +, W- =- 10 + 4- V-. 

c+ + c_ C+ + C- c+ + c_ C+ + C- 

Together with the known formulas V-(t, 2 ) = 00(2 — C-t) in region (II), 10 +(t, 2 ) = 
100(2 + c+t) in region (III), we get 

o+(t,0) = 2C+ v 0 (~c-t) + C+ io 0 (c+f), 

c+ + c_ C+ + c_ 
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_ c. | — c _ 

w-(t,0) =-- 1 — W 0 (c+t) + - —Vo 

C+ + C_ C+ + c_ 

Propagating the boundary information along characteristic lines, we obtain 

v(t,x) = v+(t ——,0) for (t,x) in region (III), 
c+ 

w(t, x) = w- (t + —, 0) for (t, x ) in region (II). 
c_ 


Collecting terms, we arrive at the following solution. 

In region (I): 

v(t, x) = vq(x — C-t), w(t , x) = wq(x + C-t); 

in region (II): 


v{t, x) = vq{x — C-t), 

2 c_ ,c+ 


c+ — c_ 


w(t,x) = —^-w 0 ( — (a: + c_<)) + 

C_|_ -p C— C- ' C + + C- 


vo{-x -C-t); 


in region (III): 


v(t,x) = —^— vo( — (x - c+t )) + -- —w 0 (c + t - x), 

c+ + c_ c+ c ++ c - 

w(t, x) = wq(x + c+t); 
in region (IV): 

v(t, x) = vq(x — c+t), w{t, x) = wq(x + c+t). 


Finally, the solution u to (5.3) is obtained by integrating ^(v + w) with respect to t. 

As a first consequence of the solution formula, we see that w is not necessarily 
continuous across the line x + C-t = 0, separating regions (I) and (II), and v is not 
necessarily continuous across the line x — c+t = 0, separating regions (III) and (IV). 
This shows the singularities emanating from the origin due to lack of compatibility 
conditions of the initial data with the jump data. 

As a second consequence, we can demonstrate the existence of a reflected singu¬ 
larity by taking data 


i>o(a:) = w 0 (x) = S(x + 1), 


corresponding to delta function initial data u±(x) = 8{x + 1), ito(x) = 0. In this case, 
the transmission condition has to be understood in the sense of limits as x —>• ±0 in 
V(Q, oo) (functions of x with values in the space of distributions with respect to t). 
Similarly, the initial data are taken in the sense of C([0, oo) : V (R)). We obtain the 
solution 


v(t, x) = 8(x — c-t + 1), w{t,x) 


S(x + c-t + 1) H——- S(—x 

v J C+ +c- K 


C-t + 1 ) 


for x < 0, containing the reflected singularity, and 


v(t, x) 


———15( — (x — c+t) + l), w(t,x) = 0 
. + c_ c+ 
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Fig. 10. Propagation of delta function at x = — 1. 


for x > 0, containing the transmitted singularity. The geometry is depicted in Fig¬ 
ure 10 (with c_ = 1, c + = 2). 

By inspection, one verifies that v + w and (w — v )/c are continuous maps of x 
with values in 0, oo), as required. Further, we find that 


1 /■* 

i(t,x) = — / (v(s,x)+w(s,x))ds 

2 Jo 

—— (h(—x + C-t — 1 )H(x + C-t + 1) + — - H(x + C-t — 1 

2c_ V c_(_ -I - C— 


2c_ 

for x < 0. On the other hand, 


1 5(—(x — c+t ) + l) = — 5{x — c+t H——) 


—< 

c_ 


C±s 
C- ' 


so that 


i(t,x) 


i 2c + 


2c_ c+ + c_ 


H ( — x + c+t ——) 
v c _' 


for x > 0. In the sector above t = 1, u(t, x) has the common value e (c++c ) 011 e if^ ier 
side of x = 0. 

It remains to be shown that the Colombeau solution U G ^([0, oo) xR) with initial 
data Uq = 0, U\ = t(6(- + 1)) is associated with the distributional solution u(t,x) 
computed above. In fact, the problem can be reduced to (5.1) with a(x) = l/c(x) 2 , 
b(x) = 1 and initial data zo(x) = 0, zi(x) = (x + l)H(x+ l)a_ with a_ = 1/cL, more 
precisely, its Colombeau version 


Ad 2 Z~d 2 x Z = 0 in f7([0, oo) x K), 

Z\ t =o = 0, d t Z\ t ~ o = l(z\) in t?(K.) 


with A = i(a). Indeed, differentiating this equation twice with respect to t shows that 
U = d 2 Z satisfies the same equation, but with initial data 

U\t= 0 = dtZ\t=o = 0 , 


d t U\ t =Q = dt(d 2 Z\ t=0 ) = = t{6{■ + 1)). 

Here we have used the fact that supp t(S(- + 1)) = { — 1} and that A is constant (equal 
to a_) there. Since z±(x) = (x + 1 )H(x + l)a_ is continuous, the quoted result from 
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[18] can be applied to the classical connected solution, showing that it is the associated 
distribution corresponding to Z. Hence U = d^Z is associated with the second time 
derivative of the classical connected solution, which equals the function u(t,x) just 
computed, as can be verified by direct calculation. 
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